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Summary — The gas flow characteristics of the dise column designed by Sterpuens and Morris 
[4] are studied by means of pressure drop determinations across the column, Water is run down 
the column with a range of gases passing countercurrent through the column. The gases studied 
are hydrogen, air, vinyl chloride and dichlorodifluoromethane. 


Correlations for pressure drop across the column, with and without liquid flowing over the 
dises, and across the empty column, are presented, which enable a correlation of data for pressure 
drop across a single dise to be made. 


The pressure drop data and gas film coeflicients of mass transfer for the disc column are 
compared in order to show the relationship between mass transfer characteristics and gas flow 
characteristics. Changes in gas flow characteristics as indicated by break points in graphs 
relating pressure drop and gas velocity were found to be consistent with changes in the mass 
transfer data. Equations which enable a prediction of the location of changes in the gas flow 
characteristics are developed. 


Data are also presented for the flooding velocities of the disc column with the system air-water. 


Résumé—Au moyen de mesures de pertes de charge au travers dune colonne a disques de 
STEPHENS et Morris, lauteur étudie les caractéristiques de Fécoulement gazeux. On fait couler 
de leau dans la colonne ot des gaz divers passent a contre-courant ; hydrogéne, air, chlorure de 
vinyle, dichlorodifluoremethane. 


L’auteur présente des corrélations pour la perte de charge au travers de la colonne, avec 
ou sans écoulement de fluide sur les disques, avec la colonne vide, et en déduit une corrélation 
pour la perte de charge correspondant au disque unique. 


L’auteur compare dans le cas de la colonne a disques, les résultats sur la perte de charge et 
ceux sur le coefficient d’échange dans le gaz de fagon a établir l'interdépendance des caractéristiques 
du phénoméne d’échange et de écoulement gazeux. On trouve que des changements dans le 
caractére de écoulement qui ne manifestent pas une brisure sur les courbes : perte de charge/ 
vitesse d’écoulement, concordent avec des altérations dans le phénoméne d’échange de matiére. 
L’auteur développe des équations prédisant les altérations dans les caractéristiques de 
lécoulement. 


Données sur les vitesses d’'inondation de la colonne a disque dans le systéme air-eau. 
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INTRODUCTION 


The dise column as introduced by Sreruens and 
Morris [4] has to date been used primarily for 
the investigation of gas absorption systems in 
which the liquid film is the major resistance. To 
evaluate the liquid film coeflicients, the overall 
coeflicients of mass transfer are corrected to take 
into account the effect of gas film resistance. 
To establish the characteristics of the column 
with respect to the gas film, Srernens and 
Morris [4] applied a liquid film correction to 
overall coeflicients obtained for the absorption of 
ammonia in water from dilute mixtures with air. 
These results were then applied to other gas 
mixtures by means of a general correlation 
applicable to wetted-wall columns [1]. [2], viz.: 


U pa PD! Pass 


Rorer [3] reported data for the evaporation 
of water and carbon tetrachloride by air in a 
dise column, These data and the gas film coefli- 
cients for the absorption of ammonia in water 
from dilute mixtures with air were well correlated 
by an equation of the same form as used by 
STEPHENS and Morris. 

The use of a general correlation applicable to 
wetted wall columns for the determination of gas 
film coefficients in a dise column is not entirely 
justified on the basis of the experimental work 
reported. No data have been published for gas 
film coefficients in a disc column for any other 
gas phase except air or dilute mixtures with air 
as the diluent. 

Theoretical and experimental developments in 
the three fields of fluid friction, heat and mass 
transfer indicate that they are all closely related. 
It was therefore considered that an investigation 
of the gas flow characteristics of the disc column 
would provide data useful in the interpretation 
of gas film mass transfer coefficients. It is 
reasonable to expect comparable changes in the 
mass transfer coeflicients and the conditions of 
gas flow. Srerpnens and Morris observed a 
change in slope of the line obtained when the 
gas film coefficients for a particular liquid rate 
were plotted against relative velocity of the air 


stream. The existence of such a discontinuity 
or “ break point” in the graph is an indication 
that the gas flow characteristics of the disc 
column also undergo a change at a definite air 
velocity. 

The gas flow characteristics of the dise column 
were conveniently studied by pressure drop 
measurements across the dises of the column. By 
selecting a range of gases which had extremely 
different physical properties it was possible to 
establish relationships which showed the positions 
of discontinuities observed in the gas flow con- 
ditions and enabled the pressure drop measure- 
ments to be correlated with the velocity and 
physical properties of the gas phase. 

A knowledge of the limiting capacity is useful 
in the design of experimental work using the dise 
column. The air velocities at which the column 
ceased to function due to entrainment of the 
liquid were determined. 


Diagram of apparatus. 
H—Calibrated orifice plate 
J—Compressed gas source 
K—Needle control valve 
Thermometer 

M— Liquid flow nozzle 

Disc supporting weight 


Fig. 1. 
A—Dise column 
B— Micromanometer 
C— Liquid rotameter 
1)— Constant head tank 
E— Circulating pump 
F—Pump sump 
G— Liquid seal arrangement 
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EXPERIMENTAL Procepure 
Apparatus 
The laboratory column used was developed by 
Sreruens and Morris [4] and consisted of 27 
carbon dises, 1-45 to 1-5em diameter, 0-44 to 
0-45 cm thick, threaded edgeways on a 16 B.S.W. 
wire. The discs were maintained mutually at 
right angles by means of a plastic cement. The 
dise assembly forming the liquid flow path was 
mounted vertically inside a 2-5c¢m bore glass 
tube through which the gas was passed. The 
manometer tappings for determining the pressure 
drop across the discs were located a_ short 
distance from each end of the disc assembly. 


Table 1, Principal dimensions of disc column 
Number of dises = 27 
Dise diameter = 1475 em 
Dise thickness = 0-45 em 
Tube diameter = 25¢em 
Mean perimeter for liquid flow = 0-122 ft. 
Equivalent diameter for gas flow = 0-052 ft. 
Free space (dry) = 89% 
Absorption surface (dry) = 01605 sq. ft. 


The pressure drop readings were measured 
using a micromanometer graduated to read to 
one hundredth of a millimetre water gauge 
pressure differential. 

The general layout of the apparatus is shown 
in Fig. 1 and the principal dimensions are shown 
in Table 1. 


Procedure 


The total pressure drop across the column was 
obtained for different liquid and gas flow rates 
for a wide range of gases. The gases studied 
were hydrogen (mol. wt. = 2), air (mol. wt. = 29), 
vinyl chloride (mol. wt. = 62-5) and dichlorodi- 
fluoromethane (mol. wt. = 121). As the investi- 
gation was centred on the gas phase character- 
istics, the liquid phase flowing over the dises 
was not varied and water was used in all runs. 
Due to the construction of the column used it 
was considered necessary to determine the 
pressure drop over the column without the dises 
in position with the same range of gases. These 


readings enabled an estimate of the pressure drop 
over the dises to be made as distinct from the 
pressure drop over the complete column, including 
column walls and entrance and exit losses. 

The maximum air velocities at which the 
column could be operated without the water 
being blown from the discs and the column 
flooding, were also investigated for different 
liquid rates. The liquid rate was kept constant 
and the gas velocity slowly increased until flooding 
occurred, 


Results 

(a) Pressure drop across column 

The data for pressure drop across the whole 
column for the range of gases: hydrogen, air, 
vinyl chloride and dichlorodifluoromethane are 
shown in Table 2. The data are well correlated 
when AP, is plotted against vp®*' T°? on a 
log-log scale as shown in Fig. 2. 
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7e 
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Fig. 2. Pressure drop across the column. 
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Data for pressure drop across the column. 


| 
| 
| 
| 


Gas 


Velocity 


St. /see. 


Table 2. 
Liquid Gas Liquid 
Rate Velocity Gas Rate 
| ft. /see. Re fi2 lb. - 
hr.ft hr. ft. 
Air Water 
408 513 (16501580 
408 690 2210 02420 | 358 
408 835 2680 OB195 
408 S83 «2835 
108 3470 0-4920 368 
WS 136 4370 0-6940 368 
408 15-6 3000-8740 68 
408 16-5 5290-9625 368 
408 177 1-0970 368 
408 20-6 6620 14430 
368 
335 1-05 337 «0226 368 
335 411 0-0266 
335 142 4560-0308 166-5 
335 175 (5620-0368 166-5 
335 2-37 761 0-0492 166-5 
335 28S | 916 0-066 166-5 
| 1190 166-5 
335 490) «(1573s 3094 166-5 
335 735 2360 2440 166-5 
335 2820 0207816655 
166-5 
188 1048 166-5 
188 630 20200-1763 
188 74000-23800 188 
30200 
| 10-7 | 34300 399688 
188 13-0 4175 5280s 88 
188 | 46200-6425 
188 | 162 | 5190 | 07500 | 
| 
Dichlorodifluoromethane - Water 
2355 B96 | 7700-1970) 
255 341 | 6630 01620 131 
255 276 «58700-1150 131 
255 216 4190 O-0818 131 
255 167 3244 00504 131 
2355 134 2608 004882 131 
255 096 1865 00328 
| | | 355 
1665 3897 | 7720 01663 | 355 
1665 334 6500 01393 355 
1665 278 5400 O-1044 355 
166-5 226 864390 00779 355 
166-5 82 | 8540 00574 855 
166-5 182-2570 00410355 
166-5 092 1780 06-0256 
Vinyl Chloride — Water 
255 $27 | 5080 | 0-1457 | 353 
255 352 4180 355 
255 303 3600 0-0862 355 
255 236 ©6280 0-0594 | 355 
255 19 2310 0-0498 355 
255 1-50 17850-0369 355 
2355 | 1605 | 0-0269 
| 


Gas 
Re 


1590 
1990 
2380 
2730 
3630 
4120 
4330 
4720 
5530 
6100 
6330 
6780 


7680 
6320 
5000 
3920 
2840 
1690 


7640 
6550 
5480 
3940 
2860 
1963 
1293 


0-0246 
0-0308 
0-0369 
00470 
0-0881 
0-1294 
0-01890 
0-02650 
OO415 
1-1320 
1-2475 
1-3225 | 


01640 i 
01232 
0-0882 
0-0613 
0-0410 
0-0239 


O23 
O17T25 
01330 
0-0882 
O-0575 
0-0369 


0-1250 
0-1006 

00-0718 | 
O-O512 
0-0410 
0-0328 j 


iquid | Gas 
Rate Velocity Gas 
Ib./- fit. /see. Re 
hr ft. 
Air — Water 

| 
231 
281 655 2100 
281 820 2630 
281 10-8 3470 
281 13-0 4175 
231 5130 
2831187 6000 
281 20-7 6650 
255 1-05 337 
255 1-35 433 
255 1-80 578 
255 1-90 611 
255 2-40 772 
255 | 2316 1015 
2355 +460 1477 
255 640 2055 
255 710 «©2280 
255 810 2760 
378 1-05 337 
378 1-35 433 
378 | 546 
378 | 745 
1085 
378 1381 
378 1860 
378 
378 | 850 2730 

_ Hydrogen Water 
255 136 630 
255 152 704 
255 | 795 
255 20-5 | 950 
255 | 23-5 1090 
255 | 265 1230 
255 29:5 1365 
255 62 | 750 
255 
235 380 1760 
255 
255 365 169 
255 | 1 | 190 
255 48 | 222 
255 73 | 338 
255 9-3 430 
255 12-5 | 579 
255 100 463 
255 14-0 | 648 


O-1187 
0-2082 
O-2870 
O-4372 
0-5930 
0-8322 
1-1100 
13350 


0-0238 
0-0266 
0-0328 
0-0390 
0-0471 
O-0737 
0-1187 
0-1826 
O-2177 
0-2784 


0-0188 
06-0266 
00850 
09-0492 
OTR 

0-1169 
1725 

02210 
02996 


| 00964 
0-1087 
0-1312 
O-1745 
0-2112 

0-3060 

0-1250 

0-3510 

0-4325 

0°3670 

0-0226 

0-0246 

0-0308 

0-0452 

0-0595 

0-0881 

00696 

01024 


VOL. 
3 


1954 


| i 
| 
| 
lb. ft.2 lb. /ft.* 
r Water 
00-1438 
) | 1926 
) 0-2566 
03080 
) 04955 
0-6060 
06600 
0-7675 
1-0030 
1-1690 
| 1-2730 
296 0-0144 
0-0239 
| | | 
S77 
693 
867 | 
1300 
» | 1797 
2230 
| 2900 
5840 
| 6810 
| 7060 | 
| 
| | 
| 
 0-0226 
5020 01636 
| 4690 | 0-1475 
4180 
8570 
2280 
1930 
1537 | 
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The three straight lines of Fig. 2 and the 
location of the “ break points” are better illus- 


trated by plotting the dimensional group a! 
p 


against gas velocity. A comparison between the 
two methods of plotting is shown in Fig. 3, in 


which AP, and = are plotted against gas 


velocity for data of one of the air-water runs shown 


in Table 2. 
1 +——420 
10 10 
| 
N 
° 
——42 
3 5 10 20 30 
Fig. 3. Comparison of methods of plotting data. 
P, 
Versus 


4P, versus v 
Air-water system ( 408 Ib. /(hr.)(ft.). 


The method of plotting = was adopted for 
vp 


all runs to determine the “ break points” but it 
was considered more fundamental to present final 
empirical correlations in terms of the pressure 
drop AP. 


(b) Pressure drop across column with zero liquid 
rate 


Flow through the dry column is equivalent to flow 
through porous media with a single fluid phase. 
Such data are best correlated by plotting a 
friction factor against Reynolds Number on 


log-log paper. 


However, it is unnecessary to use the complete 
Fanning friction factor for purposes of illustrating 
the flow characteristics of a particular column. 


The dimensional group nd has been plotted 


against Re in Fig. 4 for data on the pressure 
drop across the dry column as shown in Table 3. 


i | | © VINYL CHLORIDE 
+ DICHLORODIFLUOROMETHANE 


a | 
300 6500 1000 2000 $000. 10000 


Fig. 4. Pressure drop across column at zero liquid rate. 


(c) Pressure drop across the empty column 


To evaluate the pressure drop due to the disc 
assembly it was considered necessary to determine 
the pressure drop due to the glass column itself 
without the dises in place. 

Correlation of the data shown in Table 4 is 
obtained by plotting the dimensional group 


= against Re on a log-log graph as shown 
p 


in Fig. 5. 
HYDROGEN 
+++ 
© AIR 
4 >) © VINYL CHLORIDE 
| + DICHLORODIF LUOROMETHANE 
+ 


mo 4 | | 
+ 
—- 
it i | 
1000 2000 3000 10000 
Re 
Fig. 5. Pressure drop across the empty column (discs 
removed). 
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Table 3. Pressure drop across the column at zero liquid rate. 


Gas | Gas | | Gas | | 
Velocity Gas AP, | Velocity Gas “Velocity | AP, 
fi. /sec. Re /ft? | fi. /sec. Re ib. /ft? ft. /see. Re | Ib. /ft.2 
Air Vinyl Chloride Dichlorodifluoromethane 
| 
1-35 417 HIT | 4960 01067 4-08 | 7850 0-1822 
1-67 337 00205 4420 00839 3-46 6730 O-1413 
4:30 1380 00800 320 8800 00635 309 6000 0-1169 
5-75 1845 0-1290 296 3520 00553 270 | 5250 0-0922 
7-2 2310 O1886 | 2510 00307 215 4175 0-0594 
217 00286 «1-28 146300 2020 O-O185 
1-42 457 00154 O68 950 0-0051 O75 1460 0-0103 
1-67 537 O-O184 35 4160 oos22 402 (7820 0-1887 
6-20 1990 0-1683 30 | 8370 00656 373 7250 0-1620 
8-20 2630 02420 2468 31900 005138 0-1313 
10-7 3430 03852 2660 00369 207 5770 0-1067 
142 4560 0-5980 1-66 1973 00226 | 253 4920 0-0800 VOL. 
16-3 5230 O-7155 1-35 1605 4160 0-0615 3 
19-2 6160 09665 0-91 1082 0-0082 1-82 3540 0-0452 1954 
20-7 6650 1-0975 +0 4760 0-1064 1-42 2760 0-0287 
22-2 7130 1-2310 3-45 4110 00821 O80 1566 0-0123 
240 7700 13975 3407 3650 00685 O57 «1108 O-0072 
2980 00430 | 
2-10 2500 0-0307 
1-78 2120 00225 
175000 
1-44 1713 
O-0185 


156 1855 


(b) Between = l2and¢d = 4 
AP, = 2-94 


(c) Below d 4 
AP, = 45 » 


(d) Corrected pressure drop due to dises 


The pressure drop across the column due to the 
dises was obtained by substracting the pressure 
drop due to the empty column (as determined 
from Fig. 5) from the pressure drop data due to 
the assembled column. 

An example of the method of calculation used 
to approximate the pressure drop across a single 
dise in a dise column is shown in Table 5. 

The estimated pressure drop AP, due to a 
single dise is well correlated by plotting AP, 


10+ 


(e) Flooding characteristics 

The “ flooding velocity” of the column with 
water flowing over the dises at various liquid 
rates and air passing up the column was deter- 
mined visually. The first visual appearance of 


against vp®*' IT’ on a log-log graph. 
The dimensional equations for this correlation 
are 


Let vp! — ¢ 


(a) Above dé = 12 
AP, = 98 > 


105 


flooding occurred when small drops of water were 
blown off the periphery of the wetted discs. The 
column became unstable at this point and as the 
gas velocity was increased or the liquid rate was 
increased, slugs of liquid were blown off the dises 
and surged upward through the column. 

Due to the unstable nature of the column when 
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Table 4. 


Pressure drop across the empty column (discs removed). 


i] 
| ow 
Velocity Gas 4P, | Feloci Gas 4P, Velocity Gas 4P, 
ft. /see. Re bb. ft? | ft. /sec Re lb. /ft.2 fi./sec. Re lb. 
Hydrogen Air Air 
8-7 02 00164 | 610 1958 00451 196 6280 0-2890 
11:3 52z 00205 2199 0-0553 23-5 7340 0-4020 
158 730 0-0369 7-60 24400 0-0676 545 0-0328 
90 880 00451 = B30 2664 0-0758 6-40 2060 0-0431 
23-0 1065 0-0595 1-73 555 0-0082 740 2380 0-0513 
232 | 1165 00675 2-50 8033 0-0133 9-10 2920 0-0717 
30-3 1400 00-0863 2-74 880 0-0154 115 «3680 0-1190 
37-0 1710 01150 1140 00225 132 4230 0-1558 
53200225 1268 00266 145 4650 0-1865 
12:3 570 00266 | 470 1509 0-0327 186 5960 0-2682 
16-3 755 | 00389 1878 0-0429 223 7150 0-3668 
17-2 797 0-0410 6-95 2230 0-0532 
19-7 913 00451 770 2472 0-0655 Vinyl Chloride 
23-0 10650-0573 8-20 2632 0-0820 | 
25-0 1157 0-067 | 2-50 803 0-0164 1-62 1928 0-0061 
30-3 1400 0-0862 3-02 969 0-0184 2-30 2737 0-014 
36-5 1690 01128 1236 0-0225 2-90 3451 0-0164 
13-8 638 0-0287 4-60 1477 0-0308 3-35 3987 0-0210 
15-6 722 0-0348 5-10 1637 0-0369 +10 4879 00307 
18-2 843 00431 35 2038 0-0472 400 0-0266 
14-7 680 00318 74 2536 0-0697 146 0-0061 
17-2 796 0-0410 8-2 2632 0-08000 184 2190 0-0103 
19-0 880 0-0451 59 1895 0-0410 2-25 2678 0-0123 
= 81 2600 0-0614 2-70 3213 0-0164 
Air 10-6 3400 0-09-44 3-12 3713 0-0185 
| 15-2 4870 01907 362 4311 0-0226 
1-77 368 O-0082 17:7 3670 02560 390 4641 0-0297 
268 860 0-014 416 4950 0-0317 
31 995 00-0184 445 5296 0-0389 
4-04 | 1297 0-0267 440 | 5236 0-0389 
5-10 1637 0-0369 | 195 | 2321 0-0113 
150 
| | 
Dichlorodifluoromethane 
| 2-75 5350 0-0246 
0-98 1907 00061 
150 2920 0-0103 
4080 0-0184 
2-93 5700 0-0328 


the flooding point is approached it was found 
difficult to determine accurately the flooding 


point. 


Data based on the visual determination of the 
flooding velocity for different liquid rates are 


shown in Table 6. 
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Table 5. Determination of pressure drop across a disc. 


Liquid rate Gas velocity Gas 4P, 4P,. ft? AP, 4P, 
Ib. (hr. ft.) fit. /sec. R, (from Fig. 5) eft? 
$95 1590 0/1438 0-0328 O-1110 0-00411 
368 6-20 1990 0-1926 0-0450 O-1476 0-00547 
368 740 2380 02566 | 00597 0-1969 0-00729 
368 soo 2730 0-3080 O-0755 0-2325 0-00853 
368 113 3630 04955 0-1200 0-3755 0-01389 
ON 128 4120 00-6060 0/1482 04578 0-01695 
368 13-5 4330 06-6600 0-1606 O-4994 0-01851 
147 4720 O-7675 0/1857 O-5818 0-02153 
17-2 5530 1-0030 0-2400 0-7630 0-02820 
19-0 6100 1-1690 2846 0-03272 
19-7 6330 12730 03030 0-9700 0-03590 
21 6780 14390 0-3380 11010 0-04075 
Table 6. Air velocities which initiate flooding. 200 i TT 
(Flooding visually determined.) meee 
& HYDROGEN 
Wb. ft./ sec. 80 + DICHLORODIFLUOROMETHANE 
60 
101 320 T T T 
141 270 T TT 
188 26-5 40 
281 22-5 | | 
268 215 
408 20-0 
| 
'9 
The correlation for pressure drop across a dise, 
in which AP, is plotted against | 
on a log scale shows that the graph obtained is y Gait 
discontinuous and is characterised by two break ia 
points over the range studied. These break points 
are at vp®*' [°?2— 4 and The | | 
positions of the break points for a particular | | 
system are not fixed exactly by these equations, ' | LJ 
as they tend to give high values of the critical | ? : 5 “0 70. 30 


velocities for the air runs and low values for the 
hydrogen and dichlorodifluoromethane runs. 
However, they do give an indication of the gas 
phase velocities, near which changes in the gas 
flow conditions can be expected. 

The existence of these critical points is better 


Fig. 6. 


illustrated — by 


plotting 2 


Pressure drop across a disc, 


af, against 
p 


gas 


Reynolds Number Re for individual runs with 


s4 
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liquid rates as parameters. Such a series of 


graphs is presented in Fig. 7. 


| 
| | 
2 +—++ TT + 
j 
At | 
| 10} > — 
+ TTT > 
-+ —+— ++ 4444 -- + tt ++ -+4 
Mw -+ + + -8 
| 
wo $00 1000 $0090 10000 
Fig. 7. 


AA dichloroditluoromethane-water = 255 1b. (hr.)(ft.) 
BB dichlorodifluoromethane-water I = 131 Ib. (hr.)(ft.) 
CC air-water = 408 Ib. (hr.)(ft.) 
DD. air-water = 255 Ib. (hr. (ft.) 
EE hydrogen-water = 255 Ib. /(hr.)(ft.) 
FF air and dichlorodifluoromethane [= (hr.)(ft.) 


Referring to Fig. 7 the lines 44 and BB are 
typical of a group of lines obtained for a particular 
gas phase with different liquid rates showing the 
existence of a break point at a Reynolds Number 
of approximately 4000, This break point corres- 
ponds to the lower break point of the correlation 
shown in Fig. 2. Lines 44, DD and EE are for 
different gas phases with the same liquid rate. 
The break points for these lines are approximately 
Re = 4000, Re = 950 and Re = 600 respectively. 
These break points all correspond to the lower 
break points of the correlation shown in Fig. 2. A 
break point corresponding to the higher break 
point of the correlation shown in Fig. 2 is indi- 
cated by the line CC at Re = 3000, 

The inclusion in Fig. 7 of the data for the dry 
column shows that the flow conditions can be 


well correlated by graphing a ‘ against Re as 
vp 


shown by line FF. However, once the column 
is operating with liquid flowing over the discs 
the gas flow characteristics at a constant liquid 
rate can no longer be correlated by a Reynolds 
Number alone. 

This can be clearly seen from the lines EE, DD 


and AA, which are for constant liquid rate with 
different gas phases. Hence it appears that the 
use of a gas Reynolds Number together with some 
functions of liquid rate as adopted by SrTernEeNs 
and Morris [4] and Roper [3] for a correlation 
of gas film mass transfer coeflicients, is not 
suflicient for a general correlation applicable to 
all gaseous systems. 

As can be seen from Fig. 7, the line for zero 
liquid rate is found to intersect the line BB for 
dichloroditluoromethane at a high Reynolds 
Number and then continues above the low liquid 
rate line. A similar effect was noticed if the low 
liquid rates of the air runs were plotted. The 
higher pressure drop of the dry column at the 
high Reynolds Number could be attributed to 
the increase in impact pressure losses on the 
square edges of the dry dises. The effect of liquid 
flowing is to round off the dises and thereby 
reduce the impact pressure losses. However, the 
liquid flowing over the discs will greatly increase 
skin friction and hence the general trend that 
pressure drop is increased with increased liquid 
rate. 

The fact that break points are observed indicate 
that flow conditions in the gas phase change at 
definite critical points. Hence it is reasonable to 
expect some comparable change in the gas film 
coeflicient under similar conditions. The data 
of Rorer [3] and Srepuens and Morris [4] 
were critically examined for discontinuity in the 
gas film coeflicient at gas velocities corresponding 
to the break points indicated by pressure drop 
measurements, 

The bulk of the data reported by Roper [3] on 
the evaporation of water in a disc column was 
for air velocities which are below the first break 
point on the pressure drop relations. However, 
data for liquid rates of 121 Ib. (hr.)(ft.) and 
81 lb. (hr.)(ft.) extend on both side of the first 
break point observed in the pressure drop relation 
investigated. These data have been plotted in 
Fig. 8 in which the pressure drop data for a low 
liquid rate of the air-water system are also shown. 

The break point in the gas film mass transfer 
coetlicient graph compares well with the observed 
break point in the pressure drop data. The data 
of Sreruens and Morris [4] for the gas film 
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Fig. 8. Comparison of pressure and mass transfer charac- 
teristics. 


@ = 121). /(hr.)(ft.) Data of Roper 

© 8b. /(hr(ft.) | For gas film mass transfer co- 
efficients in the evaporation of 
water in streams of air. 

+ I = 166-5 1b. /(hr.)(ft.) Pressure drop data for air- 
water system. 


coetlicient in the absorption of dilute ammonia-air 
mixtures by water, when plotted in the same 
manner, also show a break point which is con- 
sistent with that observed for pressure drop data. 


ACKNOWLEDGMENTS 


Acknowledgment is made to Professor J. P. 
Baxter and Mr. G. H. Rover, School of Chemical 


Engineering, N.S.W. University of Technology, 
for helpful discussion and interest in this work. 
Acknowledgment is also made to the Monsanto 
Chemical Company ( Aust.) for the provision of the 
research scholarship which made this investigation 


possible. 
NOTATION 


D = equivalent diameter for gas flow, ft. 
D,, = diffusivity, ft.? /hr. 
ky = gas film mass transfer coefficient, Ib. /(hr.)(ft.*) 
(atm.) 
M, = molecular weight of the diffusing gas 
P py = logarithmic mean of the partial pressure of the 
inert gas at the phase boundary and in the bulk 
of the gas stream, atm. 
P = total pressure, atm. 
4P = pressure drop, Ib. /ft.* 
4P, = pressure drop across the whole column, Ib. ft. 
4P, = pressure drop across the empty column (dises re- 
moved), Ib. /ft.? 
4P, = pressure drop due to a single dise in a dise column, 
Ib. /ft.? 
Re = Reynolds Number for the gas phase, Dup/ 
v = velocity of the gas stream, ft./sec. 
I’ = liquor rate, Ib. /(hr.)(ft.) 
8B = a function of liquid rate 
» = absolute viscosity of the gas stream, lb, mass. /(ft.) 
(hr.) 
p = density of the gas stream, Ib. /ft.* 
pq = density of the diffusing gas, Ib. /ft.* 
= the dimensional function 
n = exponent of Reynolds Number 
m = exponent of the Schmidt Number 
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Hydroextraction IX: Studies of flow through incompressible cakes 
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Summary —Hydroextraction and filtration tests on rigid beds of particles, consisting of either 
sintered metal or ceramic material, have proved conclusively the mechanism and flow equations 
for hydroextraction postulated in previous communications. 


Résumé—Les auteurs ont procédé a des essais de centrifugation et de filtration sur les couches 
rigides de particules constituées soit de métaux frittés soit de matériaux céramiques. Ces essais ont 
définitivement prouvé lexactitude du mécanisme et des équations imaginés dans leurs communica- 


tions précédentes sur la centrifugation. 


Previous work has shown that the flow mechanism 
of hydroextraction is covered adequately by 
equations proposed from this laboratory [1, 2]. 
The discrepancies between hydroextraction per- 
meabilities derived from these equations and 
filtration permeabilities measured without dis- 
turbance of the centrifuge cake are still significant, 
though probably due to causes discussed in the 
previous communications. To confirm the validity 
of the hydroextraction equations, by attaining a 
close agreement with filtration relations, tests 
were justified on cakes which were incompressible 
and invariable in permeability under the experi- 
mental range of centrifugal fields. The following 
account describes tests on porous ceramic cylinders 
which represented complete cakes and on sintered 
metal plugs representing a radial element of a 
cake. 

Tests have been made using metal dust, e.g. 
zine, to make incompressible beds. The charac- 
teristics of such beds are satisfactory in purely 
centrifuging studies, but they are difficult to test 
in filtration within the hydroextractor basket. 
Such beds were thus of limited value in the 
present work. 


CeRAMIC CYLINDERS 


Experimental Technique 
The ceramic cylinders were of a form directly 


comparable to the hydroextractor cakes dis- 
cussed earlier. The hydroextraction tests for 
permeability AK, were obtained by carrying the 
cylinder between top and bottom circular plates 
which provided the appropriate ends of a hydro- 
extractor * basket’ (Fig. la). The cylinder was 
spun about a vertical axis by a direct-coupled 
motor, the speed of rotation being measured by 
a D.C.-generator tachometer on the motor shaft. 
The control and measurement of the apparatus 
were capable of the same accuracy as those on 
the centrifuges used in the associated tests [2]. 
Water was fed to the spinning cylinder as an 
instantaneous charge of specified volume V,, 
and the time ¢’ measured for the vertical liquid 
surface to attain the inner surface of the ceramic 
cylinder as the water drained through the 
cylinder wall. 

The analytical technique involving the measure- 
ment of drainage time has been described as 
Method III earlier [2]. It is useful so long as 
two major requirements are met. The addition 
of the initial charge V;, at t = 0 must be rapid 
but not slow the basket appreciably. The liquid 
will soon attain the selected radial velocity of 
the basket but the time required represents an 
error in the following analysis and it must be 
small compared to the total drainage time, t’. The 
time t’ must be measured accurately and detection 
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A 
T 
Dal 
Fig. 1b. 
Fig. 1. Ceramic cylinder: (a) in position for hydro- 


extraction test; (b) in position for filtration test. 
X, height of cylinder ; r,, outside radius ; r,. inside radius ; 
B, cylinder ; C, one of four connecting bolts ; D, driving 
shaft; F, felt packing; G, vent; H, feed tube of cross 
section a = 0 754.m* ; ho, height of water in 1 above top 
surface of cylinder. 


of the coincidence of r, with r, when the liquid 
surface attains the inner ceramic surface is best 
carried out with a photoelectric detector [5]. The 


shielded light ray from the lamp L in Fig. 2 is 
reflected away from the shielded photoelectric 
cell C, when a water layer is maintained at 
r, <7, When the ceramic surface is not 
covered by water, it reflects as a diffuse reflector 
and the maximum illumination reaches the cell. 


Fig. 2. Elevation and plan showing lamp L and photo- 
electric cell C in position for detection of critical wetting 
condition when ry = 


The graphs in Fig. 3 exemplify the data obtained 
from a sensitive galvanometer taking the barrier 
cell output, the ordinates being in arbitrary 
units. As the water layer becomes thinner as the 
liquid charge drains through the ceramic wall, the 
cell output increases until a sharp breakpoint is 
attained. “ Rounding-off” at the break occurs 
only when the surface is covered by gelatinous 
deposit, as mentioned earlier [5]. The cell 
output may rise slightly after the breakpoint and 
this may be attributed to further air-drying of the 
ceramic surface, after the liquid layer at r, has 
entered the ceramic (r, > 7). The photoelectric 
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method gives a clear breakpoint when r, = r, 
and is far more reliable than visual observation, 
particularly in tests requiring immediate recog- 
nition of the critical wetting condition. The 
galvanometer readings were taken at appro- 
priately spaced time intervals after the addition 
of the liquid V,, and the time ¢’ for attaining 
the critical condition taken from the extrapolation 
of the lines around the breakpoint. 


DETECTOR SCALE READING 


n 
TIME FROM ADOING Vi min 

Fig. 3. Typical result of a photo-electric detection 
experiment. Break points are the critical times. 


The analysis of the experimental data for 
various volumes V, and rotation speeds ” was 
based on the flow rate equation 

loge To, Te 
This equation is valid so long as a liquid layer 
remains on the cake, Le. r, <r, and the liquid 
fills the pores of the cake throughout the cake 
thickness (r, —r,). The liquid is in streamline 
motion and viscous resistance to this motion is 
the only significant resistance. 

Equation (1) obtains at any moment from 
t=0tot=—t, assuming that throughout that 
time the liquid layer rotates at the basket speed 
n, i.e. does not take an appreciable time to be 
accelerated after the addition of V, and rotates 
at the same angular velocity at all points in the 


(1) 
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fluid. The assumption is made that the pores of 
the cake are already filled with liquid and do not 
require filling with liquid from V,. Thus for 
accuracy when the volume of the cake is com- 
parable to the volume V,, the technique should 
include the attainment of the critical wetting 
condition (r, = just prior to the addition 
of the charge V;. 

The momentary drainage rate at time ¢ is 
related to V, the volume of liquid between radii 
r, and r,: 


K, 
dt pg. log. re 


(V.+V,) (2) 


=y¥(V.+ V,) (2a) 
where y = log. 
= w X (r,? — 12) 
whence : 
t 
dt 
log, “| = (3) 


For a given speed of revolution » and cylinder 
dimensions, the terms y and V, are constant. 
Thus a plot of log(V, + V,) against ¢’ should 
give a straight line of gradient y and intercept 
of log V,, when t’ — 0 corresponding to the use 
of a zero charge volume, V;, = 0. The gradient 
y 2-303 may be analysed to find A,. For various 
rotational speeds, y should be proportional to 
n®. The line obtained from a plot of log y against 
log n, thus should have a gradient of 2-0 and an 
ordinate when n = 1 of log (47° K_X log, r,/7,). 
This plot, or its equivalent, can thus be used to 
obtain an average K, from a series of sets of 
tests at different n, each set being at constant n 
and having variable V;,. 

To test filtration permeability K,. the ceramic 
cylinders were held between top and bottom 
plates with appropriate connections for water 
input and air venting to make sure that an air 
lock did not form in the head of the space within 
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the cylinder (Fig. 1b). The techniques used 
included the measurement of flow rate to maintain 
a constant head of liquid above the cylinder, and 
also the measurement of the rate of fall of water 
in a vertical feed tube above the cylinder. 


WW FAG’ 


Fig. 4. Diagrams illustrating the derivation of an expres- 
sion for filtration through a cylinder (equations 4-9). 


The conditions at a depth h’ below the free 
surface of the water in the feed tube are shown 
in Fig. 4. The flow through the cylinder is assumed 
to fill all the pores of a material of constant 
permeability and atmospheric pressure is assumed 
at the vertical wall on the outside of the cylinder. 
Flow through the annulus of the cylinder is 
assumed to be radial and horizontal. This is a 
reasonable approximation in the present tests, 
where with the heads fh, of 35-150cem and 
approximate cylinder dimensions of (7, — 7.) > 
lem, then dp dh is only 1-3°, of Jp dv. Fora 
precise solution the flow through the cylinder 
walls would be best analysed by relaxation 
methods. The desirable accuracy of 1°, in the 
present investigations does not warrant such 
endeavour, though it is proper to focus attention 
on data obtained at high heads, where it is most 
reasonable to neglect )p M in comparison with 


dp or. 


The radial flow through the annulus of thickness 
w at radius r is assumed to be given by an 
equation for streamline motion 


¢=- Ky dx (4) 
wr 
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The flow resistance due to frictional force is 
assumed to be the only significant resistance. 

The head h available to produce flow varies 
from h =h' atr=r,toh=Oatr=—r,. Then 
for the annulus at constant depth h', with ¢ 
constant at all radii 


dr Be "0/Te (5) 
whence 
2a K,.h’. da 
6) 


Then to obtain the total flow rate q from the 
whole cylinder over its height X, considering 
br as 5h’, we may integrate equation (6) between 
the limits h’ = h, and h’ = (hg + X): whence, 


xX 
NK X 
= d = F om 
0 
(4 + (7a) 
where 


The flow rate from the cylinder thus varies 
linearly with the mean head (hy + X/2), and 
from an arithmetic plot of q against (hy + X/2) 
the permeability may be obtained from the 
constant gradient B. 

If the flow is allowed without a feed to maintain 
a head h,, then the equation (7) is appropriate to 
an instant as the head h, falls. If the head tube 
has a cross-sectional area a and no water is fed 
to the tube, the water level and hence hy will 
fall as time ¢ increases from ¢ = 0 at hy = h, to 
t—t' at hy = hy. 


Then : 
X 
= = (8) 
whence : 
B 
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Thus using a calibrated glass tube above the 30$5-—>— T 
cylinder, K, may be obtained from B by plotting so 
the semi-logarithmic equation (9) for measured : | 
t’ values as head h, falls to various h,, in the 3-00}—+—+ 7 
form of log (h, + X/2) against t’. As h, is the 
same for all the values of h, and then the | 
data should lie on a straight line of gradient 295|—Vj = 200ce L. j/ A 
— B/2-308 a, Wit 
EXPERIMENTAL Data // 
290 4 
The two ceramic cylinders tested [6] are described & 
in Table 1, and it will be noted that their per- |RUNN® SYMBOL 
meabilities should differ considerably according | = 
to the pore sizes. | 3 
The hydroextraction results for Cylinder I are 
shown in Fig. 5. The lines representing the data (THEORETICAL) e yo 
for various V, at constant should pass 
through the log V,) axis at ordinate BREAKPOINT TIME, min 
4 log V, which is known. The actual lines concur Fig. 5a. 
3 
290} 
90 
Z 
+ 
Y + 
270 —log Vo (THEORETICAL) 
iS 20 25 30 33 40 
BREAKPOINT TIME, ——» 
N® SYMBOL 7” 
a 1930 Fig, 5b. 
4 e 1610 Fig. 5. Data from hydroextraction tests, (a) on cylinder I, 
(b) on cylinder II. 
7 1480 


at an intercept where V_ = 658 ce instead of 
the theoretical 676 cc, and were accepted as 
satisfactory. The gradients y of these and similar 
lines are shown as a function of n on Fig. 6, 
confirming the expected power 2-0. From the 


mean line shown, the permeability K, was found 


to be 2-46 x 107 g/sec?. 


The constant head filtration data shown in 
Fig. 7 were obtained using a rotameter to measure 
the rates g to maintain a steady head (h, + X (2) 
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+ —+—__ 


° 


+ 


GRADIENT, y 


Agreement between experimental results and 
equation (2). 


Fig. 6. 


in the glass head tube, the rotameter discharging 
freely into the head tube. The gradient B of the 
arithmetic plot of q against (4, — XY 2) was 
constant, confirming the assumption of a stream- 
line flow in which the only significant resistance 
was that due to viscous resistance. From the 
gradient B in equation (7a) the filtration per- 
meability was found to be 2-16 « 107 @ sec?, 


cc/sec 
° 


~ 


° 


FLOW RATE, g 
4 


° 


° 20 40 60 80 100 120 140 cm 
(to + 


Agreement between experimental results and 
equation (7). 


Fig. 7. 


Table 1. Data on Ceramic Cylinders* 


Cylinder I 
Manufacturer's grade S.P. (Sand pot) 
Nominal maximum pore size 35 
Mean dimensions : 
height, X 
outside diameter 
inside diameter 


Cylinder 
P6A Porcelain 
2 microns 


15-20 
12-22 


15-32em 
14-98 em 
13-41. em 
7 496m 
570 ce 


4903 
Cylinder volume ; 


676 


* The ceramic cylinders used were made by Doulton and 
Co. Ltd., Albert Embankment, Lambeth, London, S.E.1. 


The filtration tests under falling head from 
h, 150 cm at 0 are plotted in Fig. 8. The 
gradient (— B 2-303 a) of the straight line from 
equation (9) gave a permeability of A, = 2-36 

107g sec*. These data from falling head 
tests are generally more reliable than those at 
constant head due to the time anomalies in flow, 
which can occur at these low permeabilities, being 
less in the falling head tests. The agreement 
between A, and K,, is probably within 5°,. 

The hydroextraction data for Cylinder IL are 
shown in Fig. 5, 6, leading to a permeability 
K, = 0-29 © 107 g sec? from the mean line on 
Fig. 6. 

The constant head filtration data for Cylinder 
II are shown in Fig. 7, the gradient of the line 
giving K, — 0-31 x 107 g sec*. 

The falling head filtration data for Cylinder I 
in Fig. 8 lead to a value of K, — 0-30 ~« 107 
g 

The data show a satisfactory agreement between 
hydroextraction permeability A, and _ filtration 


Agreement between experimental results and 
equation (9). 


Fig. 8. 
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permeability A. confirming the equations pro- 
posed. Since the ceramic cylinders readily 
available were limited in dimensions, porosity 
and evenness of porosity distribution, further 
tests were carried out with sintered metal plugs. 
It was hoped that these would provide data in 
closer agreement than those from the ceramic 
cylinders in view of the success of small cell 
experiments (3) wherein the liquid level could 
be maintained at r, — 0, independently of 
time. This technique was expected to be superior 
to the time tests used on the ceramic cylinders 
(6), and it was expected that with small metal 
samples more uniform material would be available 
in greater widths (7, — 7.) in a wide porosity 
range. It should also be emphasised that in 
tests of the present type the cylinders must be 
regarded as * expendable.” After relatively short 
time of immersion in water algae growth and 
slight gelatinous films on the ceramic made the 
cylinders unsuitable for further work. All water 
used in the batch tests was treated with either 
thymol or mereuric chloride but some of the 
continuous tests required volumes of filtered 
sterilised water beyond the laboratory capacity. 
The series of tests on small plugs were expected 
to be very much quicker and thus more satis- 
factory than the long series discussed above. 


Sinrerep Merat 


The sintered metal plugs used in these tests are 
described in Table 2. The material consisted of 
small metal spheres moulded sintered 
together without much deformation. The plugs 
were produced without machining the circular 
end faces of the rods. The cylindrical surface 
was sealed by filing and soldering when necessary. 
The bulk of the plug was not affected by this, the 
circular ends remaining fully permeable, appro- 
priate to the pore dimensions, and flow being 
possible only through the length of the plug. 
These plugs were held in an appropriate carrier, 
which is shown in Fig. 9a. The carrier was whirled 
about the vertical axis at a constant speed of 
rotation. Flow rates q were controlled to maintain 
the water level at G in the feed tube over the 
experimental range of speeds x. The carrier was 
then held with the plug vertical (Fig. 9b) and 
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Table 2. Data on Sintered Bronze Plugs 
Plug No. Length, L Diameter, d_— Sinter grade* 
cm cm 

1A 380 175 A 
1B 380 1-75 B 
1c 3-80 1-75 
2A 204 1-21 A 
2B 210 1-20 B 
2c 210 1-20 ( 
BA 5-08 0-90 A 
5-00 OSS 
4A A 


* The sinter grades are classified as follows according 
to the manufacturers, Sintered Products Ltd., Sutton Cold- 
field, Notts, England : 


Maximum particle to pass 


Grade (thousandths of an inch) 
A Ol 
B O2 
Cc OS 


the filtration test was conducted under suction 
with a constant feed of distilled water to keep the 
water level at G. 


° ° 
G 
8 
~ p 
a b 


Plug carrier, (4) in position for hydroextraction 
test ; (b) in position for filtration test. 


Fig. 9. 


A, bronze plug; B, plug jacket; the common surface 
between A and Bis sealed ; C, a perforated base cap tightly 
screwed on B; I1, rubber gasket ; D, tapered brass con- 
nection tightly screwed in B; E, rubber tube connecting 
D and glass feed tube F'; G, level of water in feed tube ; 
L, driving shaft ; K, counterbalancing weight ; 1JN, back 
brass plates connecting B, L, and K, and securing tube F 
(the similar front plates are not shown) ; M, rubber tube ; 
P, to suction pump. 
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The hydroextraction permeability was assessed 
as A. in the flow equation : 

To — Ve 


The filtration permeability A, was calculated 
from the corresponding equation : 


AH .K, 


11 


The relation between A, 
directly assessed as : 


H 


Kp q P. 


and K, could be 


(12) 


where P, = centrifugal head = 3. 


Metal plugs with permeability of the order of 
107 g sec® showed a continuous decrease in flow 
rate when tested for filtration under constant 
hydrostatic heads. The application of suction 
below the plug reduced the effect but did not 
eliminate it. The flow rates might decrease to 
one fifth of the original within a few minutes in 
spite of using distilled water and maximum 
precautions to ensure cleanliness. Use of electro- 
lyte solutions did not alter the results. Vibration 
and intermittent suction of air had a large effect, 
even up to temporary elimination of the pheno- 
menon. This vibration may be the reason why 
the time effects are far less noticeable in hydro- 
extraction, 

Plugs of permeability of 10% g sec? showed 
negligible time effects when tested under the same 
hydrostatic conditions over a considerable length 
of time, whereas with similar plugs of permeability 
of 10° g sec?, the effect was so quick that it was 
impossible to obtain reproducible results with 
aqueous solutions. With K = 107g sec? the 
test had to be done as quickly as possible to 
attain the initial maximum flow rate, whereas 
plugs of permeability 10° and 10% g sec? were 
more convenient to test. Flow rates through 
plugs of permeability 10 and 10° g sec? were 
so high that the hydroextraction tests were 
difficult due to air being sucked into the feeding 
tube. The range of K around 10° to 10° g/sec* 


was the optimum for comparison between per- 
meabilities in hydroextraction and filtration with 
the present apparatus. 


Experimental Data 

Fig. 10 is a typical plot of the q, n relationship 
for the case of plugs of permeability near 10° 
g sec, The gradient 2-0 was also obtained with 


FLOW RATE 
: 
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Fig. 10. Agreement between experimental results and 
equation (10). 
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Table 3 Relation between K, and K » for Sintered Bronze Plugs 


| Filtration Test 
Plug No. SS 
| P Kr x 108 
g/sec* 
1A | 73-0 660 2.02 
1B 73-4 25-20 7-61 
2B 2040 7-38 
2c 7146 37-00 20-7 
750 


3c | 735 
| 


‘tae K,/ Ky 
n’ K, = 108 | 
r.p.m cc/min. g/sec® 
1350 7-80 2-08 1-03 
1260 25-10 7-70 1-01 
1120 48-00 | 18-68 1-00 
1250 20-40 | 7-45 1-01 
960 33-00 20-5 0-99 
1330 8-40 12-9 1-02 
Mean 1-01 


* ry = 10-50em 


all plugs of other permeabilities given in Table 2. 
Clean connecting tubes and distilled water are 
essential in these tests, as the least contamination 
on the surface directly leads to erroneous results 
and lower q, gradients. 


c | 
—+——- + + — 
| 
—+—- --+ -- +—— ——__ + 
| | | | 
| 
a | 
= | 
| + + + 
° 10 20 30 40 SO 60 70 cmHg 
SUCTION HEAD, P ——> 
Fig. 11. Agreement between experimental results and 


equation (11). 


Fig. 11 illustrates the linear proportionality 
between q and H for a filtration test on a plug 
of permeability 2-07 « 10° g sec?. The linearity 
could not be illustrated for lower permeabilities 
because of the time effects discussed above. Thus 
the range of permeability 10° and 10° g ‘sec? 
was chosen for comparison data for performance 
under hydrostatic and centrifugal heads. Table 3 
exemplifies the results obtained. The ratio KK, 
is very near to unity with an average difference 
less than + 1%, 
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CONCLUSION 

Though the agreement between AK, and A, for 
the ceramic cylinders was only within 5%, it is 
conceivable that the discrepancies were due to 
inaccuracies in technique rather than in the 
postulated mechanism for equation (1). The 
probable lack of uniformity in the cylinders, 
assumptions of negligible acceleration times in 
the batch tests, and similar approximations led 
to the 5°, discrepancy, just as the hydroextraction 
reproducibility on Fig. 6 is shown to be + 5°,. 
The method used for the measurement of A, 
with the metal plugs is more accurate than the 
batch technique [2]. 

The agreement to 1°, between K, and A, with 
the metal plugs is considered to prove conclusively 
the mechanism of hydroextraction embodied in 
equation (1). The mean differences of 10-20%, 
between A, and K, in tests on cakes of chalk 
and starch in hydroextractors and in a small 
filtration cell [3] are thus likely to be caused 
by variations in cake porosity and _ interface 
resistance as postulated earlier [4]. 


NOTATION 
A = cross-sectional area of plugs (em?) 
a = cross-sectional area of feed tube in filtration tests 
(em?) 
B = constant in equation (7a) 
gravitational acceleration (cm, sec*) 

H = hydrostatic head in filtration tests on plugs (em 

water) 


M. M. Hanunt ef al. Hydroextraction IX 


h = hydrostatic head in filtration tests on cylinders (em 


water): A’ at r 
limits of h. 


re hg at XN = 0, hy and hy are 


= permeability ; A. from hydroextraction tests, Ap 


from filtration tests (g sec?) 

speed of rotation (rev see) (r.p.m.) 

pressure (dynes em?) 

centrifugal head (em water) 

rate of flow; q¢° flow through a circular clement of 
cylinder, q, and qp for flow under centrifugal and 
hydrostatic heads respectively (ce sec) 
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Heat and mass transfer at high humidities in a wetted wall column 
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Summary~— Simultaneous heat and mass transfer data were obtained from the operation of an 
adiabatic wetted-wall column with countercurrent flow of air and water. The runs were made 
over a limited range of Reynolds numbers 2390 to 9095 and from 0-03 to 0-85 mol fraction of the 
diffusing vapour in the gas film. The data were compared with the equation of previous investi- 
gators obtained at low humidities. Correlation of the data by means of a Nusselt-type equation 
guve 

(kg = 0-021 SOM 


at a liquor rate of 77 Ib. hr. This equation differs from the GrLLILAND-Suerwoop equation 
mainly by the exponent 0-83 on the group pyyy P. Comparison with the Cunrox-CoLsuRn 
equation showed that their heat and mass transfer factors did not apply at high humidities. 
Correlation by new heat and mass transfer factors gave 


=jp = Se? (pay / 
= 0-025 


The theory of CotnurNn and Drew for the effect of simultaneous mass transfer on the heat 
transfer coeflicients was unable to explain the data. 


Résumé— Résultats des mesures simultanées d'échanges de chaleur et de matiére dans une 
colonne adiabatique a parois mouillées avec écoulements opposés d'air et d'eau. 


Expériences avec des nombres de Reynolds variant entre 2390 et 9095, et des concen- 
trations molaires pour la vapeur diffusant au travers du film gazeux variant entre 0,03 et 0,85. 
Comparaison avec les formules proposées par d'autres expérimentateurs pour de faibles humidités. 
La corrélation conduit a formule du type de Nusselt : 


(kg RT d/D)(Ppyy/ PP — 0.021 SOM 


pour un débit liquide de 77 Ib. /hr. x ft. Cette équation différe de celle de Gritann et 
par lexposant 0,83 pour le groupe pyy,/ P. La comparaison avec léquation 
de Cniwron ~ CoLBuRN montre que leurs facteurs d’échanges pour la chaleur et pour la 
matiére ne sont plus applicables aux humidités élevées. 


La corrélation des nouveaux coefficients d’échange conduit aux égalités : 


in = (h/eG). 
= ip = ke Pam Se? (pyy/ 
0.025 


La théorie de Corsurn et Drew pour leffet dun transfert simultané de matiére sur le 
coefficient de transfert thermique ne peut expliquer les résultats des auteurs. 
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INTRODUCTION 


A limited number of papers have reported 
experimental results for mass transfer occurring 
at high concentrations of the diffusing vapour. 
The mean mol fraction of the diffusing vapour 
has not exceeded 0-5 in any of the investigations 
previously published. The first paper to examine 
the effect of these high concentrations is that of 
Hanks and McApams [9] who reported the 
measurements of Hemron [10] of the rates of 
absorption of ammonia from air by water in a 
single contact chamber stirred at 350 r.p.m,. in 
both the liquid and gas phases. A series of runs 
were made in which the total pressure of the gas 
was constant (1 atm) and the composition of the 
gas in the chamber was varied between 0-06 and 
0-65 mol fraction ammonia, Hanks and McApams 
correlated these results on the basis of the theory 
and the Stefan equation 


N, DPAp RT rp yy, (1) 


where « is the stagnant gas film thickness. The 
liquid film resistance was neglected, so that the 
gas film coeflicient was assumed to be equal to 
the overall coeflicient, viz. 


DP Ap 


Ka = ke RT op yyy 


The experimental coetlicients were correlated by 


ke 02357 Puy (3) 


The correlation is unsatisfactory, but it does 
indicate that k, is influenced by p,,, P. It is 
probable that the liquid film resistance is not 
negligible and that the reaction in the liquid 
appreciably affects the liquid film coetlicient. 
Hixcutey and Himes [11] varied the ratio 
Pam P from 0-55 to 0-99 in their study of the 
evaporation of water from open pans by natural 
and forced convection. ©The transfer 
coetlicient for foreed convection was not affected 
by variations in p,, P. The results for natural 
convection indicate the mass transfer coeflicient 
to be approximately inversely proportional to 
this ratio. Wape [15] evaporated seven different 
liquids from a pan 3-5 in, square into a tunnel 
472 in. 236in. The evaporations were made 


os 


into still air (natural convection) and into air 
streams 4-4 to 13-1 ft. per see. The ratio py, P 
was varied between 0-68 and 0-99. The gas film 
coeflicient, k, was independent of py, P for 
forced convection, and was inversely proportional 
to Pyy P for natural convection. 

Borvrer, Gorpon and Grirein [2] working 
with a 1 ft. diameter pan found that the mass 
transfer coeflicients for the evaporation of water 
into air at Latm total pressure were best corre- 
lated by 

ke = (4) 


for natural convection, They also correlated the 
data by the dimensionless equation 


k, RTd D 0-648 (5) 


in which the ratio of p,,, P does not enter. Their 
data may be better correlated by the dimension- 
less equation (6) than by equation (5), viz., 


k, RTd D — 0-61 (Gr.Se.pyy, (6) 


The ratio of p,,, P was varied from 0-5 to 0-99 
and the water temperature from 75°F to 202°F. 
The intluence of pyy, in equation (4) is due to 
the close relationship between the modified 
Grashof number, Gr, and pyy,. and is not due 
to the effect predicted by the film theory and the 
Stefan equation. 

Neither Hixcuiey and Himus, nor Wapr give 
suflicient experimental details to recalculate their 
data in terms of the dimensionless equations, so 
that it is not possible to decide whether their data 
for natural convection can be correlated in accord- 
ance with equations (5) or (6). 

GILLILAND and Suerwoop [8] studied the 
evaporation of nine different liquids in a wetted- 
wall tower 1-05 in. id. operating at pressures 
from 0-145 to 3-06atm. They correlated their 
experimental results on the basis of the Stefan 
equation (1) and the film theory. For Reynolds 
numbers greater than 2000 and a constant liquor 
rate of 790 cc min. (104 1b hr. for water) their 
results were correlated by the dimensionless 
equation 


(kgRTd/D) P) 0-028 Reo (7) 


No variation of the left-hand side of equation (7) 
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with pressure was found, from which the authors 
concluded that k,, was inversely proportional to 
Pum 28 indicated by the Stefan equation (1). 
However as Suerwoop and Picrorp [14, p. 77] 
point out, the ratio p,,,/P was varied only 
slightly. Hence no conclusions may be drawn 
for the effect of p,,, and P on ky. 
and Suerwoop also measured the heat transfer 
coeflicients for heating air in upward and down- 
ward flow in the column without liquid circulating. 
The experimental heat transfer data agreed 
closely with the McApams [12, p. 168] equation (8) 
for heating air in turbulent flow in a pipe 


Nu — 0-023 Re®*® Pro (8) 


Their mass transfer data are about 35°, greater 
than their data for heat transfer. 

Barner and Kose [1] studied heat and mass 
transfer in a wetted-wall column 1-049 in. id. For 
the evaporation of water into air streams at low 
humidities, their mass transfer data were corre- 
lated by the Sherwood-Gilliland type equation 


D)(pyy/P) 0-016 Re (9) 


at liquor rates varying between 10 and 100 Ib hr. 
Taking the value of the Schmidt number under the 
conditions of their experiments as 0-60 equation (9) 


becomes 
D)(pyy P) 0-020 Sc? (10) 


Their mass transfer values are approximately 13°, 
less than those of GiLuiLaANp and SHerwoop ; 
about 3°, out of the 13%, difference could be 
attributed to the different liquor rates used. Their 
mass and heat transfer data were also correlated 
in the form of the Chilton-Colburn dimensionless 
groupings, as follows 


G) Se? ® — 0-025 Re? (11) 
and 
(h/eG) — 0-031 Re®? (12) 


The heat transfer data without simultaneous mass 
transfer are 24°, greater than the mass transfer 
data. 

CotpurN and Drew [6] have considered the 
effect of the sensible heat carried out by the 
diffusing molecules during simultaneous heat and 
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mass transfer. Based on the film theory and in 
the absence of condensation in the gas phase, they 
developed the relation 


=h = ha (1 — e*) (13) 


where h, is the heat transfer coeflicient in the 
absence of simultaneous mass transfer. The 
dimensionless group, a, is the ratio between the 
nett sensible heat carried by the diffusing vapours, 
and the heat transfer in the absence of mass 
transfer. For one diffusing component, a_ is 
w,¢, hy and is positive if the mass transfer is in 
the same direction as the heat transfer. For two 
diffusing components, is (w,e, + hy. Uf 
the heat transfer is in the same direction as the 
mass transfer, h is greater than h, and conversely 
if the heat transfer is in the opposite direction to 
the mass transfer h is less than h,. 

There is still some confusion in the literature 
concerning the effect of the inert gas concentra- 
tion on the rate of diffusion of a vapour through 
a gas film. The effect of p,,,/P indicated by tlhe 
GILLILAND-SHERWOOD correlation and by the 
CHitron-CoLBURN correlation is retained purely 
on the basis of its appearance in the theoretical 
equation (1). There appears to be no definite 
attempt in the literature to show the influence 
of large variation of p,,, P on the rate of mass 
transfer, and in view of the importance of p,,,/P 
in such cases as the condensation of vapour from 
non-condensable gas [4], [7] and in the con- 
densation of mixed vapours, this investigation 
was made, 


Apparatus, Procepure anp 


An outline of the apparatus is given in Fig. 1. The 
column and procedure were similar to those used 
by GiLuiLanp and Suerwoop. Air from a blower 
was passed through a_ glass-wool filter, and 
metered by a Fischer and Porter rotameter. The 
humidity of the incoming air was measured by 
wet- and dry-bulb thermometers before it entered 
an electric preheater where its temperature was 
controlled by a variac. The temperature of the 
air leaving the preheater was normally in the 
range 220-240°F. After preheating, the air was 
mixed with steam from the boiler. The boiler, 
of 6 gallons capacity, was heated with a 3-kW 
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240V AH RW RECORDING WATIMETER 
nov 
Fig. 1. Outline of the experimental apparatus. 


Helicoil A.G.E. twin element immersion heater 
sealed into the side of the boiler and was lagged 
with lin. asbestos sheet. Power input was 
controlled by a variae on each element and 
recorded on an Evershed and Vignoles recording 
wattmeter. The steaming rate of the boiler was 
determined from the power input and a cali- 
bration curve determined under the conditions 
of the experiments. A correction was made for 
the condensation in the line from the boiler to 
the point where the steam and air mixed. The 
steam-air mixture was further heated electrically 
to the temperature at which it entered the column. 
The wetted-wall column was a glass tube 0-90 in. 
id. 87} in. long. The column was enclosed in 
another glass tube 1-3 in. id. around which was 
wound an clectrical heater. The column heater 
was controlled by a variac adjusted to maintain 
the operation as adiabatic as possible as deter- 


mined by a heat balance. The heater was enclosed 
by a third glass tube. ‘This arrangement allowed 
the water film to be seen at all times so that a 
check was kept to ensure that complete wetting 
oeceurred during all runs. The temperature of 
the water leaving the separator at the base of 
the column was measured. The water from the 
separator passed to a vessel with a controlled 
electrical heater. The water level in this vessel 
was maintained constant by additions from a 
100 ml. burette. The water was recirculated by a 
small centrifugal pump, the water temperature 
entering the top of the column being adjusted 
to the temperature at the base by the controlled 
heater. The rate of water recirculation was 
measured by a rotameter. 

The air-vapour mixture from the top of the 
column passed to a condenser except for those 
runs in which the load on the condenser was 
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excessive. The product which was condensed was 
collected and measured. The humidity of the 
cooled air-vapour mixture was determined by wet- 
and dry-bulb thermometers. The temperatures 
of the air-vapour mixtures entering and leaving 
the column were measured by copper-constantan 
thermocouples and a Cambridge thermocouple 
potentiometer. The lines at these two points 
were well lagged so that the radiation correction 
was negligible. 

Manometers at the inlet and outlet of the 
column were used to determine the total pressure 
in the column. 


Mernops oF CALCULATION 
The runs were interpreted by comparing the 
experimental values of dimensionless groups with 


the values predicted by the correlations of 
GiLLILAND and Suerwoop, Cron and 
and and Drew. The experi- 
mental data is shown in Table 1. The water 
vapour entering the column is the sum of the 
water vapour from the boiler and in the incoming 
air corrected for the condensate from the steam 
line. The water vapour leaving the column in 
runs 1-19 was calculated from the amount entering 
the column and the amount added to the cir- 
culating water (make up). In runs 20-27 the 
water vapour leaving the column was calculated 
from the humidity of the exit air and the con- 
densed product. In these runs the amount of 
make-up was small and the laticr method was 
considered more accurate, The rate of evaporation 
was determined from the make-up in runs 1-19 


Table 1, Experimental data 


Dry air Boiler Con- Product Make Mixture Temp, Column Liquor Water Temp. 
Run Air lb. lb. dry air rate densate rate uprate at Column Pressure rate at Column °F 
No, rate Ibs. hr. lbs. hr. lbs. hr. lbs. hr. atm, Ibs. — 

‘Ibs. hr. In Out In Out In Out 

1 (WTS O1S 2-39 Os 2-76 O51 629 336 1-007 77 183-9 1841 
2 OOT 2-39 2038 O78 Oly B82 1-008 83 | 1743 | 1744 
3 1734 O07 156 77 175-4 175-6 
4 O10 +73 O12 140 O34 383 77 1820 | 182-8 
5 775 OO 573 Ol OAT 1-06 370 77 188°8 | 189-0 
6 “ooo “O16 +73 O17 OST 630 379 1007 77 198-2 198-0 
7 2359 O12 113 61s 351 1-008 77 1641 1645 
OOS 2-39 1-41 B52 1-008 77 1589 158-7 
” 2440 2-39 O13 1-77 338 1-005 77 152°3 151-9 
10 189 “Olt O24 651 S808 442 1-008 151 204-9 204-7 
225 “25 584 O24 657 O-ST 777 435 0-987 136 202-5 202-8 
12 292 OOF 028 584 O19 6-70 112 810 0-992 111 200-8 200-9 
13 334 “O22 6-85 1-16 813 454 0-992 11 200-0 200-0 
14 479 OOF O12 678 116 766 437 0-995 105 195-7 195-6 
15 572 007 5340 677 1-28 424 1-001 ill 193-3 193°5 
16 524 OOF O12 6-87 137 1-000 195-2 195-0 
7 434 “O31 O16 O85 1-29 S806 1-000 111 197-1 197-3 
18 669 06 “O42 6-79 1-30 736 421 1-000 111 191-2 191-8 
632 006 “O89 O12 6-30 1-32 761 4350 1-000 192-5 192-4 
20 770 O15 Oo oOo O14 - B02 184 1-007 24 101-9 101-9 
21 771 oo oOo 0-16 - 301 71 1-007 143 102-0 102-0 
22 771 oo O15 - 297 176 1-005 42 101-6 101-6 
25 772 Oo Oo O16 297 172 1005 1020 
25 673 oo oo O12 O17 | 271 156 0999  T7 95 98-3 
26 11-50) 0083 “033 Oo oo oo 263 167 77 1025 102-6 
27 2407 Ooo oo oo O34 205 145 77 
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and from the difference between the entering and 
exit water vapour in runs 20-27, The inlet and 
outlet humidities were calculated from the water 
vapour entering or leaving the column and the 
dry air rate. The driving force, Ap, was taken 
as the difference between the vapour pressure 
corresponding to the liquid temperature and the 
partial pressure of water vapour in the gas phase. 
In caleulating the gas film coeflicients of mass 
transfer, the logarithmic mean between the 
driving foree at the top and bottom of the 
column was used. The true mean driving force 
as determined by graphical integration differed 
by no more than 1°, from the logarithmic mean 
in any of the runs. The values of the gas film 
coetlicient of mass transfer were calculated from 


In evaluating the Reynolds number the arith- 
metic mean of the mass rate of flow at the top 
and bottom of the column was calculated. The 
viscosity was determined at the mean film con- 
ditions by the method of Bromiey and WiLKe [3] 
using the viscosity data of Perry [13, p. 370]. The 
diffusivities were calculated from Spaupina’s 
equation [14, p. 10] for the same conditions. The 
temperature used in the dimensionless groups 
was the mean film temperature. The value of 
Puy Was taken as the arithmetic mean between 
the logarithmic mean partial pressures of air in 
the gas film at each end of the column. The 
Schmidt numbers were calculated for each run, 
and as the variation was slight an average value 
of 0-55 was taken. In calculating the dimension- 
less groups for correlating with p,,,/ P?, the values 


ke = Ny A Apia. of these groups were corrected to a liquor rate of 
Table 2. Calculated mass transfer data 

H at column | | | 
Run Ib. Ib. dry air Na | Re Reo PM HM). Reo? Mass 
No. lb. atm. bp \\D P balance 

In Out moles hr at 77 Ibs. /hr. | at 77 Ibs. 
1 O578 ©7070-0283 0-043 | 0-90 2389 18-7 0-0294 0-469 0-0210 + 78 
2 O828 0439 06-0433 0-061 0-97 3417 204 0-0236 0-587 0-0219 — 38 
3 0339 0-428 0-0867 0-071 167 35-4 0-104 0-583 0-0187 
4 0610 0-059 181 6520 374 0-0253 0-501 0-0209 
5 O735 O87T2 O-0589 0046 183 5308 38-2 0-0310 0-413 0-0206 + 66 
6 1482 1655 06-0483 0-036 1-84 3987 37-6 0-0386 0-269 0-0166 — 92 
7 0286 00628 121 5100 255 00214 0-685 0-0228 
8 OOTRS OO8T 123 6826 26-1 O-O1T2 0-735 0-0199 
9 06-0996 0-172 O-0983 166 9095 35-3 0-0183 0-0234 
1 «2074 0-020 350 3187) 0-0760 0-158 0-0172 - 
1) «2493 2880 0-025 264 S273) «499 0-0555 O-175 0-0146 | —172 
12) 08-0623 3511) 0-0651 0-212 0-0216 + 
13 20630-0644 0-034 3662) 0-0511 0-230 0-0185 = 26 
16 1200 0644 0-0438 0-301 0-0209 34 
15 1007) «1231 O-OTI2 239 4531 468 0-0413 0-337 0-0222 + 68 
16 1361) 07620-0404 0-0455 0-0228 + 58 
17 1316 «1613 «O-OTIT 0-:0357 275 3995 52-6 0-0515 0-280 0-0228 + 16 
18 1057 06-0723 252 4828 45-4 0-0380 O-371 0-0226 0-0 
19 O91) 1120 00734 0-0481 209 4641 406 0-0350 0-355 0-0197 + 28 
20 «06-0068 60-0338 06-0323 0-496 2814 126 0-0187 0-951 0-0273 
21 «40-0075 06-0363 ©0128 06-0297 0-500 2822 15-0 0-0187 0-949 0-0271 
22 08-0350 OO117 0-554 13-6 0-0196 0-950 | - 
28 0-0363 O-0122 00297 0-562 2831 144 O-O18T 0-0270 
24 «608-0077 06-0349) 06-0811) 06-516 2831 131 0-0179 0-950 0-0260 - 
25 0-0312 06-0095 00290 0-446 2508 11-6 0-0175 0-955 0-0254 
26 06-0083) «06-0330 «60-0345 «(0-619 16-1 | 0-0156 0-950 0-0187 + 
27 06-0062 06-0212 06-0200 06-0237 116 8925 313 0-0165 0-966 0-0244 59 
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77 Ib. hr. from the relations of Fig. 2. The mass 
balances are expressed as an error in closure 
defined by 


®, error in closure 
_ total vapour in — total vapour out 


100 
make-up 


for those runs in which the complete data could 
be obtained. In some runs the condenser was 
unable to handle the load, and in others the 
amount of make-up was too small to measure 
accurately. 

The calculated mass transfer data are given in 
Table 2. 

Only those runs in which the error in closure of 
the heat balance was less than 10°, were used in 
calculating the heat transfer data. For these runs 
the heat transferred across the gas-liquid inter- 


face is the product of the mass of water evaporated 
and its latent heat of vaporisation at the water 
temperature. Some of this heat (never more than 
15°.) is transferred by radiation and the rest is 
transferred by convection. The heat transferred 
by convection was taken as the difference between 
the total heat transferred and the heat transferred 
by radiation as calculated by the approximate 
method of Perry [13, p. 490). In calculating 
the convective heat transfer coeflicients, h, the 
log-mean temperature difference was used. The 
true mean temperature difference by graphical in- 
tegration differed by less than 2°, from the log- 
mean in any of the runs, The specific heat data for 
airand water vapour were taken from WENNER [16]. 
The variation of the Prandtl number was slight 
and the average value of 0-76 was taken in the 
calculations, The dimensionless group, a, derived 


Table 3. Calculated heat transfer data 


| | | | | | 
Run qT Stat St. Pri, Re? Heat 
No. BTU hr. hr. F. | h TT Ibs. /hr. at 77 Ibs. | hy | balance %, 

2 | 75 737 277 3-65 W535 W227 0-91 0-126 — 35 

1589 1512 279 741 0227 0-91 0-124 | — +4 

4 1383 1336 310 3-89 0213 0-85 0-132 —73 

5 1044 991 481 0203 0-81 | —43 

6 | 852 788 $73 | 0-76 O115 | —17 

- | - - | | | - 

9 1782 | 1764 205 S17 0249 100 O40 — 53 
10 - - - | - 
- 7 ~ - | - - 
12 | 1094 973 400 333 | 00381) 0163 + 43 
13 1134 1010 408 339 0160 064 0148 + 22 
14 1137 1039 383 S72 066 O186 +382 
15 1257 £1168 | 366 4°36 00413 0186 O74 56 
16 | 1344 1241 304 431 00418 0186 0-74 0-155 —70 
17 1264 1152 w2 3892 O71 | | —O7 
18 1278 1189 366002 W399 0-72 0-139 0-0 
19 1297 1202 38000 00395 O178 0-142 0-0 
20 27 217 | 182 225 00555 227 O91 0-055 0-0 
21 238 | 123 26500556 0227 0-91 0051 +47 
22 217 | #125 | 287 0568 W232 0-93 0-053 + 52 
23 28 0-050 +17 
24 217 124 240 00346 0223 O89 0-050 + 58 
25 176 0 0-046 + 58 
260 200 105 B78 00576 0256 102 0049 | 82 
375 | 81 633 00469 O22 | 0-035 8-4 
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by and Drew is equal to — 18N Ahy 
for the uni-directional diffusion taking place in 
the column. The heat balances were calculated 
for a datum of 100°F and are based on 1 Ib. of 
dry air, The balances are expressed as an error 
in closure defined by 


®. error in closure 
(heat in — heat out) 100 
latent heat of water evaporated 


A sample calculation (run 13) is as follows : 
Dry air rate — 3-34 Ib. per hr. 


Total vapour entering column 
0-08 — O14 
5-73 lb. per hr. 


Humidity at column inlet 
5-73 3-34 — 1-716 th, per lb. 


Rate of vaporisation 
make-up rate — 1-16 Tb, per hr. 
0-06.44 Ib. moles per hr. 


Humidity at column exit 
(5-73 1:16) B34 
2-063 Ib. per Ib. 


Vapour pressure of water at mean temp, (200-0°F) 
0-783 atm. 


Total pressure in column 
0-992 atm. 


Partial pressure of water in gas phase 
at column inlet 0-729 atm. 
at column outlet 0-763 atm. 


Log mean driving force 
0-084 atm. 


Mean driving force by graphical integration 
0-034 atm. 


Inside area of column (dry) 
0-731 sq. ft. 


Mass transfer coetlicient ky, 
00-0644 0-034 0-731 
2-59 Ib. moles (hr.) (ft.?) (atm.) 


at column inlet 
(0-992 — 0-729) (O-092 
In (0-263 0-209) 
= 0-236 atm. 


— 0-783) 


Pray ®t column exit 
(0-992 — 0-763) — (0-902 — 0-783) 
In (0-229 0-209) 
0-219 atm. 


Mean 
(0-236 4 0-219) 2 


0-228 
Mean 
P0992 


Mean mass rate of flow W 
4+ 3-7? 116 2 
9-65 Ib. per hr. 


Film temp. 
at column inlet 
at column outlet 
mean value 


(818 4+ 200) 2 — 507 F 
(454 + 200) 2 — 827 F 
(507 + 827) 2 417° F 


Viscosity at mean film temp. 
water vapour — 00165 centipoise 
air 0-0252 centipotse 
~ mixture 0-0185 centipoise 
0-0448 Ib. (ft.) (hr.) 
by the method of Bromiey and Winker 


Cross-sectional area of column (dry) 
= 0-00442 sq. ft. 


Mean Reynolds number for gas phase flow 


9.65 3-142 0-0888 (0-90 12) 
- 
0-0448 


D = 1-46 10% P(T + 441) 
0-992 (417 + 460 + 441) 
= 2-55 sq. ft. per hr. 


k,RTd/D 
2:59 O-T2Z8 (417 + 460) (0-90 12) 
48-6 


Liquor rate 
= 111 Tb. per hr. 


Correction factor (from Fig. 2) 
1-048 


104 


R. C, Catens and G. H. Rorrer : Heat and mass transfer at high humidities in a wetted wall column 


(k, RT ( Re®**) at 77 Ib. per hr. 
48-6 1-048 908 
0-0511 
Mean molecular weight of gas stream 
20-7 
Se calculated for this run 


Se used in calculations 
0-55; Sc?* — 0-671 
(keP yy M,, Se? at a liquor rate of 77 Ibs. hr. 
2-59 O228 20-7 0-00442 
9-65 1-048 


0-671 
= 0-00358 
(ke Pum M,, G) (Se?) (Re®?) = 
= 000858 5-16 — 0-0185 


Mass Balance 
Total water vapour in 
as humidity in air (3-34 . 0-008) 
0.03 Ib, per hr. 


from boiler (5-84 — 0-14) 5-70 ,, 
make-up 1-16 ,. 
total Oso 


Total water vapour out 
condensed product 6-85 Ib. per hr. 
in outlet air (3-34 « 0-022) 0-07 ,, 
total 6-92 ,, 


Error in closure 
(6-89 — 6-92) 
= — 26% 


100 


Total rate of heat transfer across gas-liquid inter- 
face 

1:16 977-7 — 1184 Btu per hr. 
Heat transferred by radiation, using the approxi- 
mate method of Perry [13, p. 490] 

= 124 Btu per hr. 


Heat transferred by convection 
= 11384 — 124 
== 1010 Btu per hr. 


Log-mean temperature difference 


408°F 
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Mean temperature difference by graphical inte- 
gration — 400°F, 


h convection heat transfer coeflicient 
1010 0-731 408 — 3-39 Btu (hr)(sq.ft.) 


¢e, at column inlet 
O246 + 1716 
~ 1-075 Btu (F) (tb. of dry air) 


ec, at column exit 
0-243 4+ O475 2-063 
1-223 Btu (-F) (lb. of dry air) 


mean value of ¢, 
1-149 Btu (Ib. of dry air) 


heG 
1-149 3-34 


Liquor rate correction factor (from Fig. 2) 
1045 


h cG at liquor rate of 77 Ib. hr. 
0-00390 1-045 


(h eG) (Pr? *) (Re) 
0-00373 0-833 5-16 = 0-0160 


hy 0-025 
a = —18N 


18 = 0-0644 « 0-479 


= — 0-148 
0-731 3-39 (0-025 
Heat Balance— 
based on 100°F datum, 1 Ib. of dry air 
Heat in 
— sensible heat of air-vapour stream 
(713 x 1-075) 766 Btu 
—- latent heat of water vapour 
(1-716 1037) 1779 Btu 
—— sensible heat in water evaporated 
(1-16 = 100) 3-34 35 Btu 
total 2580 Btu 


(813 — 200) — (454 — 200) 
“In (613/254) 
|| 


Heat out 
sensible heat of air-vapour stream 
(354 1-223) 
latent heat of water vapour 
(2-063 10387) 2139 Btu 
total 2572 Btu 


433 Btu 


Error in closure 
» 
- 2580 2572 100 
2139 — 1779 


22%, 


Discussion 


The effect of varying liquor rate on the mass and 
heat transfer coeflicients was determined from 
several runs operated at low humidities and at 
almost constant Reynolds numbers for the gas 
stream. The results of these runs are shown in 
Fig. 2 where the dimensionless groups, h eG and 


+4 YT 0080 
> 
oon 
| 
2 46 66 66 100 20 140 160 
WATER RATE Ib/ hr 
Fig. 2. ‘The influence of the water rate on the mass and 


heat transfer coeflicients. 


k,RTd DRe®™ are graphed against liquor rate. 
The liquor rate has a small but appreciable effect 
on the gas film coeflicients of heat and mass 
transfer. 

The results were plotted by several of the 
methods suggested by previous investigators. The 
values of dimensionless groups involving &,, or h 
were corrected by means of Fig. 2 to a standard 
liquor rate of 77 1b. hr. and Suer- 
woop [8] correlated their results by the Nusselt- 
type equation (7). For a Sehmidt number of 0-55 
and a liquor rate of 77 Ib. hr. the equation may 


be reduced to 


(kA, RTd D)(pyy P) 0-017 (14) 
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The experimental mass transfer coeflicients at low 
humidities agree with the GILLILAND-SHERWoOoD 
equation, but, as shown in Fig 3, the data at 


Comparison of the mass transfer data with the 
GILLILAND-SHERWOOD equation (14). 


Fig. 3. 


high humidities fall consistently below this line. 
The data at low and high humidities from this 
investigation were correlated by the method of 
least squares to give 


(k,RTd D)(pyy, — 0-016 (15) 


or 
(k,RTd D)(pyy PY — 0-021 Se (16) 


at a liquor rate of 77 lb. hr. and with a standard 
deviation of 10°,. Comparison of equations (14) 
and (15) show that the present data do not sup- 
port the usual procedure of assuming that k, 
is inversely proportional to p,,, P. The data were 
also correlated using viscosities and diffusivities 
evaluated at the bulk conditions of the air- 
vapour stream. The equation derived is almost 
identical to equation (15). 

and [5] have correlated heat 
transfer data in an empty tube by plotting a 
heat transfer factor, (h/cG)Pr?*, against the 
Reynolds number. Similarly, they correlated 
mass transfer data by plotting a mass transfer 


— | 
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factor, (kgPyM,,/G@) Sc?*, against the Reynolds 
number, Their correlations may be written 
(h PP? (ke pyyM,,/@) 

0-023 (17) 
For a liquor rate of 77 lb./hr. equation (17) be- 
comes 


(h/eG) = (keP yyy M,,/G@) 


0-026 (18) 
010 
— 
4 


CHILTON ~COBURN 


8 
eltee 


Comparison of the mass transfer data with the 
equation (18), 


Fig. 4. 


The mass transfer data of this investigation have 
been correlated in Fig. 4 by plotting (kgp..,,M,, G) 
Sc?® Re? against P. The Chilton-Colburn 
equation (18) is represented by the horizontal 
line. The data fall consistently below this line 
at high humidities although the data at low 
humidities agree with it. The data are best 
correlated by defining a new mass transfer factor, 


Jp = @) (pyy/ PY? (19) 


Jp is found to be equal to 0-025 Re-®? at 77 Ibs. /hr. 
liquor rate. Comparison of equations (18) and (19) 
show that the present data do not support the 
Curtron-CoLBuRN mass transfer factor at high 
humidities. 

The heat transfer data have been correlated in 
Fig. 5 by plotting (h cG) Pr?* Re? against 


010) 
2 CHILTON - COLBURN 
Nip 
INVESTIGATION 
0-01 | 
O10 100 
—— 
Fig. 5. Comparison of the heat transfer data with the 
equation (18). 
Pay P. The Cuiron-CoLpurn equation (18) is 


represented by the horizontal line. At low 
humidities the data agree with the Cuiron- 
CoLBURN equation. At high humidities the data 
fall consistently below the line of equation (18). 
The data are best represented by defining a new 
heat transfer factor, 


Ju = (h eG) PY??? 


Jn is found to be equal to 0-025 Re®? at 77 
Ibs. hr. 

Both the heat and mass transfer data of the 
present investigation can be correlated by the 
one equation 


in = 0-025 Re®? 


at 77 Ib. 

The correlations above apply to all the data 
irrespective of whether the heat balance is 
positive or negative. Therefore any small tem- 
perature gradients in the liquid due to slight heat 
leakages are negligible. 

CotpurN and Drew [6] have proposed a 
theoretical equation for heat transfer coeflicients 
with simultaneous mass transfer. The effect of 
the simultaneous mass transfer on the heat 
transfer coefficient is given by equation (13), 
where a, in this case, is — 18N,¢,/Ahy. The 


(20) 


(21) 
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equation (13) of Connurn and Drew. 


CONCLUSIONS 


ties. 


ratio h h, is plotted in Fig. 6 as a function of a. 
The correlation is poor. 
and run 9 show that the correlation in Fig. 5 
against Pp», P is superior. The theoretical line 
of equation (13) is shown in Fig. 6, The predicted 
effects of a are much smaller than those obtained 
experimentally, 


In particular, runs 2-6 


+ + | + 
| 
+ 
| 
| 
| | 
i 
Fig. 6. Comparison of the data with the theoretical 


The liquor rate in a wetted-wall column has an 
appreciable effect on the mass and heat transfer 
coeflicients in the column. 

The gas film coetlicients of mass transfer in a 
wetted-wall column vary inversely as the 0-83 
power of the ratio p,,, P. and not as the 1-00 
power as has previously been assumed on the 
basis of the film theory. 
and Suerwoop, and of and 
CoLBuRN therefore should not be used at high 


The equations of 


The gas film coeflicients of heat transfer in an 
adiabatic wetted-wall column vary as the 0-2 
power of the ratio p,,, P and are not independent 
of this ratio as required by the Cuiurox-CoLBuRN 
equation. 

The theoretical equation of CoLaurn and Drew 
based on the film theory is unable to explain the 
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reduction in the heat transfer coetlicient, occurring 


at high humidities, found in this investigation, 
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NOTATION 


a — the ratio between the nett sensible heat carried 


by the diffusing vapours and the heat trans- 
ferred in the absence of mass transfer 

the area of the gas-liquid interface, taken as the 
inside area of the column, when dry, sq. ft. 
specific heat, (Ib.) ( F.) 

humid heat of the air-vapour mixture, Btu. (Ib.) 
of dry air) ( F) 

specific heat of water vapour, Btu (Ib.) ( F) 
the specific heat of diffusing component 1 or 2 
respectively, Btu (Ib.) ( F) 

diameter of column, ft. 

molecular diffusivity of the diffusing compo- 
nent, sq.f:. hr. 

base of Naperian logarithms 

gravitational constant ft. (hr.*) 


mass velocity, Ib. (hr.) (ft.?) 


modified Grashof number for mass transfer 
(gd® p* (M, T, M; — 1) 
convective heat transfer coeflicient, Btu, (hr.) 
(ft.*) (F) 
convective heat transfer coeflicient in the ab- 
sence of mass transfer, Btu (hr.) (ft.2) (/F) 
humidity, lb. of water vapour per Ib, of dry air 
mass transfer factor — Se? 
(Pum pyeu 
heat transferfactor =(h eG) Pr? (pay 
thermal conductivity, Btu (hr.) (ft.) 


- the gas film coeflicient of mass transfer, Ib. 


moles (hr.) (ft.2) (atm.) 

the overall mass transfer coefficient, Ib. moles 
(hr.) (ft.2) (atm.) 

molecular weight at the gas-liquid interface 

the molecular weight of the gas-vapour mix- 

ture “ far away “ from the gas-liquid interface 

the molecular weight of the gas-vapour mixture 

in the main stream 


- the rate of evaporation of water vapour, Ib. 


moles, hr. 
Nusselt number — hd k 
total pressure, atm. 


- the logarithmic mean partial pressure of the 


non-diffusing gas in the gas film, atm. 
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dp — the driving force for the diffusing vapour, atm. ; = absolute temperature at the gas-liquid inter- 
4p, the logarithmic mean between the driving forces face, “R 
at each end of the column, atm. T, = absolute temperature “far away” from the 
Pr Prandtl number = eu/k gas-liquid interface. 
a the rate of heat transfer by convection, Btu br. 4t = temperature difference, °F 
qr — the total rate of heat transfer, Btu /hr. (41). ~ logarithmic mean temperature difference, F 


W), Wy — the rate of diffusion of component 1 or 2 res- 
pectively, lb. (hr.) (ft.*) 
W — mass rate of flow, Ib. /hr. 


gas constant, (ft.*) (atm.) (Ib. mole) 
Re Reynolds number dG, 


Se Sehmidt number — «/pD w — distance in the direction of diffusion, ft. 
St Stanton number = A /cG viscosity, Ib. (hr.) (ft.) 


T absolute temperature, R p density, Ib, /ft.” 
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Variables affecting activity of molybdena-alumina hydroforming 
catalyst in aromatization of cyclohexane 


C. G. and Warson 


University of Wisconsin, Madison, Wisconsin 
17 February W54) 


Summary— An exploratory study has been made to determine the nature and magnitude of 
the more obvious variations in the activity of a single batch of catalyst. A pelleted, copre- 
cipitated, gel-type catalyst simular to that used in many commercial catalytic reforming units 
was studied in a differential reactor. 


Apparently the molybdenum oxide promoter was at a higher level of activity when the oxide 
was in a reduced state, the improvement being directly proportional to the degree of molar 
reduction of the oxide. Deposition of coke on the catalyst, which reduced activity exponentially, 
appeared to be a surface-controlled reaction and to involve a dual site mechanism. Increased 
hydrogen to hydrocarbon molar feed ratios reduced carbon deposition in direct proportion to 
the associated increase in hydrogen partial pressure but also reduced conversion because of the 


lowered reactant partial pressure. 


1954 


Catalyst regeneration appeared to occur in two distinct steps. In the first step, the rate had 
a second order dependency on coke concentration. In the second step, for which first-order 
dependencies on coke concentration and oxygen partial pressure were observed, mass transfer 
seemed to be the controlling mechanism. No improvement in catalyst activity was noted in 
the removal of carbon in this second step. 


Résumé—Avant dentreprendre étude du meécanisme réactionnel, on a procédé a étude 
préliminaire de la nature et de la grandeur des variations les plus importantes observables avec une 
méme préparation de catalyseur. 


Dans un convertisseur différentiel, on a étudié un catalyseur semblable a celui qui est utilisé 
dans de nombreuses unités de “ reforming “ : catalyseur pastillé & partir d'un gel obtenu par 
coprecipitation. 


Il semble que lactivité promotrice de Foxyde de molybdéne est plus levee quand loxyde est 
dans son état réeduit, Fameélioration paraissant directement proportionnelle au degré de réduction 
molaire de Toxyde. Le dépét du coke sur le catalyseur réduisant lactivité exponentiellement 
semble étre une réaction contrélée par la surface avee un meécanisme faisant appel a deux sites. 
Avec des alimentations correspondant a un rapport molaire aceru de Phydrogéne a Vhydrocarbone, 
le dépét de carbone est réduit proportionnellement 4 Taccroissement de la pression partielle 
de Thydrogéne, et également 4 la diminution de la conversion qui est associ¢e & Tabaissement 
de la pression partielle de lhydrocarbone. 

La régénération du catalyseur parait s‘effectuer en deux ¢tapes: Dans la premiére, la 
vitesse dépend de la concentration du coke suivant lorde deux. Dans la deuxiéme étape, of 
lon observe des relations de premier ordre avec la concentration en coke et la pression partielle 
de Toxygéne, il semble que le transfert de matiére est le mécanisme contréleur. L’enlévement 
du coke dans cette deuxiéme étape ne s'accompagne pas d'ameélioration dans lactivité catalytique. 


* Present address : Explosives Department, Eastern Laboratory, E.1. du Pont de Nemours & Co., Gibbstown, 
New Jerev, U.S.A. 
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The investigation was divided into the following 
studies : 


(1) Reduction of catalyst promoter by pre- 
treatment with hydrogen for different lengths of 
time and in various states of dilution with nitro- 
gen, with and without water vapour addition. 

(2) Catalyst fouling by coke deposition during 
hydrocarbon processing — its effects on catalyst 
activity, a suggested mechanism of deposition, 
and its suppression by hydrogen. 

(3) Catalyst regeneration — burning rates, 
degree of regeneration, and mechanism of burning 
during the later burning stage. 


Marertais 
Catalyst 


coprecipitated molybdena-alumina gel-type 
hydroforming catalyst in the form of jy in. eylin- 
drical pellets was used. It was obtained as a sample 
from California Research Corporation, identified 
by them as CRC-7242, and was similar in chemical 
composition to that used in many of the catalytic 
reforming units in the petroleum industry. It 
was specified as containing approximately 10°, 
molybdena as MoO,, as having a specitic surface 
in the order of 90m? gm, and as being well 
below 1°, in impurities, chiefly iron and sodium. 


Feed 

A 98 — mole per cent cyclohexane, purchased 
from Shell Chemical Co., was distilled in a 30-plate 
Oldershaw column at a reflux ratio of 7 1, the 
top 5°, and the bottom 5°, being discarded. 
Negligible discrepancy was observed on an 
analytical infra-red differential spectrum between 
this feed stock and a sample of very carefully 
purified cyclohexane. 


Reactor 


The most significant feature of the equipment was 
the differential reactor used, that is, a reactor 
containing only a small catalyst bed which pro- 
vided conversions generally no greater than 10°,,. 
Thus, the temperature and partial pressure 
variations were most probably either negligible 
or sufliciently small over the catalyst bed to 
permit the use of linear approximations. 
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The 30 grams of catalyst were contained in a 
lin. stainless steel pipe equipped with a pre- 
heater spiral to vaporize cyclohexane and to 
promote mixing and heat transfer. A well, con- 
taining thermocouples at levels near the top, at 
the middle, and near the bottom of the 3 in. high 
catalyst bed passed axially through the centre 
of the bed. Because of the large thermal inertia 
of the aluminium-bronze blocks in which the 
reactor pipe was contained, the temperature of 
the catalyst bed could generally be maintained 
within 1 F. 


Feed and Recovery Systems 


Cyclohexane was pumped into the reactor by a 
calibrated assembly consisting of a Zenith gear 
pump which discharged a heavy oil. This oil 
displaced mercury which in turn displaced cyclo- 
hexane from a high-pressure burette. Copper 
tubing was used for liquid feed lines. Hydrogen 
feed was controlled by a needle valve and 
measured by a high-pressure orifice. In order to 
eliminate the oxygen contamination, all hydrogen 
was passed through a stainless steel tube con- 
taining quartz chips heated to 1100 F, and then 
through Drierite, before passing to the reactor. 

In the recovery system, ice water was circulated 
through the liquid product condensers, U-tubes 
containing Drierite and Ascarite, in that order, 
were used for water and carbon dioxide deter- 
mination in product gases. In order to permit 
complete quantitative determination of carbon, 
the regeneration gases were passed through a 
tube of CuO particles held at 300°C, to oxidize 
any carbon monoxide, 


EXPERIMENTAL Procepure 

Catalyst Selection 

Catalyst pellets from the original sample lot 
described were selected for the one experimental 
catalyst bed used throughout the investigation 
on the basis of their crushing strength. This was 
considered necessary in the light of previous 
experiences in which the greater fraction of the 
pellets in the catalyst bed were observed to soften 
and crumble. Furthermore, it was also considered 
desirable because such a selection tended to fix 
an important variable in the catalyst effectiveness 
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factor. Thus only pellets were used which with- 
stood a 15-lb. crushing force applied perpendi- 
cularly to the cylinder axis. This selection test 
was passed by only 22% of the pellets in the 
batch of catalyst sample originally received and 
which was presumably representative of that 
commercially used. It has been indicated that 
it is probably reasonable to assume direct 
proportionality between crushing strength and 
density. The pellets in the batch finally selected 
had densities between 1:47 and 1-70, with an 
average of 1:56 gm ce, compared to a density 
range of between 1:13 and 1-70 for the original 
batch, with an average density of 1:36 gm ce. 


Reaction Conditions 
Studies were made at the following conditions : 


Variable Range “Standard” 


Catalyst weight, W (gm) 312 B12 
Reciprocal space velocity, 

WF (gm catalyst gm 

mole C,H,, hr.) 5-22 10-4 
Pressure (atmabs.) 1-7-7-8 17 
‘Temperature 860-980 O40 
H, C,H, molarfeed ratio, 1-71 


Run Procedure 


The “ standard ~ run consisted of the following 
operations performed in this order : 


(1) Catalyst pretreatment with 100°, Hy, at 
0-09 sefm. and 1000°F, 15 min. 
(2) Hydrocarbon processing, 15 min. 


(3) Purge with N, at 0-25 sefm., 5 min. 

(4) Regeneration with air at 1-7 atmabs. and 
at 0-025 sefm.. 15 min. 

(5) Purge with N, at 0-25 sefm., 5 min. 


The liquid product sample used for the deter- 
mination of catalyst activity was that obtained 
during the 3 to 10 min. interval of hydrocarbon 
processing. As will be noted from Fig. 3, the 
decline in catalyst activity during this period 
was still rapid. However, the rate of decline was 
considerably lower than during the first 3 min., 
so that reproducibility in sampling was con- 
siderably enhanced. The samples were analyzed 
for benzene concentration by a checked infrared 
spectrophotometric method, 


Determination of Catalyst Activity 

Catalyst activity factors were determined by 
reference to standard curves of conversion, 2, vs. 
reciprocal space velocity, WF, with tempera- 
ture as a parameter, as suggested by Houcen 
and Warson [4]. These curves, illustrated in 
Fig. 1. were obtained in a single series of runs 
using the same batch of catalyst which was used 
throughout the investigation. The “standard” 
procedure and conditions (except temperature, 
which was a parameter) were used. Since the 
standard procedure included sampling during 
the 3 to 10 min. hydrocarbon processing period, 
the catalyst activity which provided a basis of 
comparison was by consequence of this procedure 
integrated over this same time period, 


Since 

L, — activity of catalyst under test condi- 
tions ; 


L, — activity of catalyst under standard 
conditions ; activity considered unity; 


(W F), — standard reciprocal space velocity at 
x, T, P (from standard curves) ; 

(WF), — experimental reciprocal space velocity 
ata, T, P; 

av, T, P — conversion, temperature, and pressure, 


respectively ; 


= 
(WY FP) 

The difference between uncorrected plots of 
conversion vs. reciprocal space velocity, and plots 
corrected to a standard carbon concentration is 
also illustrated in Fig. 1. These plots, as they 
have been used in this study and as generally 
used for the determination of catalyst activity, 
apply to a standard process period. Hence, each 
point represents a different integrated catalyst 
activity because of the varying coke concentration 
on the catalyst resulting from the variations in 
operating temperature and in the conversion 
corresponding to the different feed rates. However 
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Fig. 1. Conversion vs. W/F: contrast between uncorrected curves and curves 
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at the lower conversions and reciprocal space 
velocities, nearly linear relationships generally 
apply. Thus the error is usually not serious in 
this range. It is obvious from the curves shown 
in Fig. 1 that a larger discrepancy is introduced 
in directly comparing catalyst activity at different 
temperature levels using the uncorrected curves. 
This results from the fact that uncorrected iso- 
therms refer to the catalyst at substantially 
different levels of activity because of the appre- 
ciably different carbon concentrations at the 
different temperatures. 

It is also desirable to observe the correction 
to a common coke concentration in obtaining 
reaction velocity constants and enthalpies of 
activation. It is often preferable to obtain the 
more meaningful values referring to the clean 
catalyst from a plot corrected to zero coke 
concentration. The conversions for the clean 
catalyst at a given temperature would be obtained 
by extrapolation. 


CaTaLyst PRETREATMENT 
A striking increase in the activity of the catalyst 
was observed to result from a hydrogen pre- 
treatment, and this was accompanied by the 
evolution of appreciable quantities of water. It 


has been tentatively suggested that the molyb- 
denum promoter was reduced to a lower and 
more catalytically active state of oxidation than 
that in which it was left after the catalyst 
regeneration with air. 

The following reaction seemed to be a plausible 
description [6] : 


5MoO, + 8H, —-» Mo,O, + 8H,0. 


The quantity of water evolved in pretreatment 
was in good agreement with that calculated by 
assuming the molybdena to be present entirely 
as MoO, after regeneration and to be completely 
reduced to Mo,O, during an extended hydrogen 
pretreatment. 

The data showed a 10°, increase in activity 
with a fivefold increase in hydrogen flow rate 
during pretreatment, from 0-018 scfm. to 0-090 
scfm. At the lower of these flow rates, the 
reduction in hydrogen partial pressure over the 
bed at the point of maximum rate of water 
evolution was practically negligible. 

The effect of the length of the pretreatment 
on catalyst activity is shown in Fig. 2. The 
greater percentage of the improvement in activity 
occurred during the early period of pretreatment. 
The fact that the curve plotting the degree of 
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Fig. 2. Effects of pretreatment time on activity and on degree of catalyst reduction. 


reduction (defined as the weight per cent oxygen 
removed, based on fresh catalyst) was of the 
same shape suggested the linear relationship 
between activity and the degree of reduction, 
shown graphically in Fig. 3. 

Changes in the degree of reduction were also 
accomplished by varying the concentration of 
water vapour in the pretreatment gas between 
0 and 2 mole per cent while holding pretreatment 
time constant. While the rate of reduction was 
considerably reduced by the added water vapour, 
the activity was again observed to be a linear 
function of the degree of reduction, as also shown 
in Fig. 3. Possibly water was adsorbed by active 
sites and thus tended to diminish the accessibility 
of the active sites for reduction. In hydrocarbon 
processing, water may then have been rapidly 
displaced from the active sites by the preferentially 
adsorbed hydrocarbons so that the changes in 
activity were a reflection only of the changes in 
the degrees of reduction. Thus, Fig. 3 also shows 
that the relationships between activity and 
degree of reduction effected by both methods 
coincide. Chemisorption of water at these 
temperatures has been evidenced in the dehydro- 
genation of butane on chromia-alumina catalysts ; 
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Fig. 3. Activity factor vs. degree of catalyst reduction. 
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and present data on the molybdena-alumina 
catalyst used in this work indicates that the 
surface of the lower oxide is saturated with 
chemisorbed water [9]. 

The effect of the concentration of hydrogen in 
the pretreatment gas diluted with nitrogen was 
also observed. The activity increased only a few 
per cent for an increase in hydrogen concentration 
from 25 to 100 mole per cent. At 0% hydrogen, 


duced during the hydrogen pretreatment. The 
reaction under test was the aromatization of a 
paraffin. Using a chromia-alumina catalyst in 
the dehydrogenation of butane, Stevens [8] 
observed an appreciable decline in catalyst 
activity with the reduction of the chromium oxide 
promoter and observed a corresponding increase 
in catalyst activity with oxidation of the chro- 
mium oxide. He was able to relate the effects 
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Fig. 4. Conversion vs. process time at various conditions of total pressure 
and H,/CgH,, molar feed ratio (all other conditions standard”). 


the effect was the same as that of no pretreat- 
ment. No data were obtained at concentrations 
between 0% and 25% hydrogen. 

The reduction of molybdenum oxides has been 
mentioned as taking place during hydrocarbon 
processing in commercial hydroforming for pro- 
duction of high-octane motor fuel [3]. The final 
metal oxide at the completion of the 4} hr. 
hydrocarbon processing period has been suggested 
to be approximately Mo,O,. This reference states 
that the reduction of the molybdena results in a 
decline (not an increase) in catalyst activity. 
and Ripeat [5] have attributed 
the promoter action of an alumina supported 
catalyst to the stabilization of the MoO, pro- 


on the catalyst activity to the partial pressures 
of hydrogen and oxygen and the time of treatment 
for the reduction-oxidation cycle. 


HyprocarBon PRocessING 

Nature of Carbon Deposition 

The initially rapid decline in conversion with 
process time is shown in Fig. 4 and is reflected 
in the plots of activity vs. process time and vs. 
carbon deposition, Fig. 5. Process time and 
carbon deposition are shown to affect activity 
exponentially, As implied here, carbon deposition 
increased exponentially with time. Within the 
experimental error, the rate of carbon deposition 
could be described as 
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dCe K 
de 


. where 


dCe 
= rate of carbon deposition, 
dé - 


Ce — concentration of carbon on the catalyst, 
and 
AK - 
Voormies [10] observed the same relationship 
in his work with both natural and synthetic 


a constant. 


was nearly ten times greater. This appeared to 
rule out the possibility that the controlling 
mechanism for this system was diffusion and 
suggested a reaction-controlled mechanism. 

A frequently advanced hypothesis of coke 
formation is the irregular polymerization of 
unsaturated by-products to form polymers which 
tend to be preferentially adsorbed and then to be 
progressively dehydrogenated. In the light of 


this hypothesis, a suggested reconciliation of the 
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Fig. 5. Decline in activity with process time and with carbon 
deposition (“ standard” conditions, except process time). 


cracking catalysts in both fixed and fluidized beds 
using a gas-oil feed. Because a_ postulated 
diffusion-controlled mechanism, wherein the coke 
represents the resistance, would imply a rate of 
carbon deposition inversely proportional to the 
coke concentration, and because of the low 
temperature coeflicient of carbon deposition which 
he observed (a 200°F increase in temperature in 
the 900°F vicinity approximately doubled the 
carbon formation rate), he suggested that a 
diffusion mechanism might be controlling the 
carbon deposition. However, the temperature 
coetlicient observed for the system used here 


high temperature coeflicient to the inverse pro- 
portionality of the rate to the concentration of 
carbon on the catalyst is indicated by the following 
analysis suggested by [9] : 
Assume that the reaction taking place in coke 
formation is 
A, + A—> Ay 
A,, ~ adsorbed coke molecule (postulated to 
be unsaturated) ; 


adsorbed 
added ; 


adsorbed larger coke molecule. 


where 


unsaturate molecule to be 


116 


Neglecting the reverse reaction and assuming that 
the mechanism of the indicated reaction involves 
reaction of an activated unsaturate on one active 
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active site, the rate of reaction is expressed as 
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at 940° F (from approximately 4°, to 16°, con- 
version) and a corresponding increase in reciprocal 
space velocity, W) F (from 5 and 22 gm gm mole 
-hr.). While a definite relationship was shown 
between aromatization 
activity, as described earlier, it is possible that 


coke deposition and 


r k coke formation, presumably a very irregular 
Raga + Ay) reaction, was promoted in part by sites which 
where did not effect aromatization. Conceivably, the 


C — a constant at the given reaction condi- 
tions ; 


a — thermodynamic activity ; 
K = adsorption equilibrium constant of sub- 
seript ; 
HC — adsorbed hydrocarbons not participating 
in reactions of coke formation, 


Assuming (1) that (A 4,a4,) >> (1 + Ayay 
K yj(4y-) on the basis that the coke is probably 
much more strongly adsorbed than the unsaturates 
or other hydrocarbons ; (2) that the coke can be 
considered in terms of thermodynamic activity 
as shown in the rate equation; and (3) that the 
activity of the unsaturate additives remains 
reasonably constant on the basis of the reasonably 
constant conversion observed during the run 
beyond the initial processing period, the rate 
equation then simplifies to 
f= 


Here 


relatively slight increase in coke deposition with 
conversion might have been the result of the 
promotion of certain side reactions producing 
unsaturates, which by their condensation con- 
tributed to coke formation. These side reactions 
may have been promoted, at least in part, by the 
site systems responsible for aromatization. The 
rate of coke deposition may also have been 
increased by the higher concentration of unsaturate 
intermediates of the principal reaction, the 
concentration of these intermediates being a 
function of the conversion. 

The hypothesis of the compound site system 
for coke formation may be further supported by 
the fact that aromatization catalyst activity was 
definitely diminished by small concentrations of 
water in the feed, but the activity for coke 
formation apparently remained relatively 
uninhibited. On the other hand, this might have 
been a result of competitive active site adsorption 
in which the coke and its intermediates competed 
against water more successfully than the cyclo- 
hexane. 


r — rate of coke deposition ~ — Again, in the aging of the catalyst, the increase 
in coke formation activity was about four times 
a4, — activity of coke ~ Ce, coke concentra- the increase in aromatization activity. On the 


tion ; 
C’ = a constant ~ A, another constant, 


and the equation is equivalent to the rate equation 
originally indicated by the data, 


dCe K 
dd Ce" 
This analysis indicates that perhaps the coke 
deposition was controlled by a dual site surface 
reaction mechanism. 
Coke deposition was observed to increase only 
about 50°,, with a fourfold increase in conversion 


other hand, this might possibly be explained in 
terms of specifically unobserved phenomena which 
occurred during the early hydrocarbon processing 
period. During this period the initial carbon 
deposition was relatively large, and yet the 
sample on which aromatization activity was 
based was not taken during this initial period but 
rather during the 3 to 10 min. period. 
HEINEMANN [2] cites evidence which might be 
interpreted as lending further support to the 
hypothesis of the compound site system in coke 
formation. In the dehydrogenation of methyl- 
evelopentane using a chromia-alumina catalyst, 
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he observed that the effect of adding a small 
quantity of benzene to the feed was to reduce 
the coke deposition to approximately a third of 
its original level without appreciably diminishing 
the conversion. This might have implied the 
poisoning of sites active in carbon formation, 
suggesting a separate site system for the dehydro- 
genation and for the coke formation. 


Hydrogen Addition 


The effect of hydrogen addition with the 
hydrocarbon feed is shown in Fig. 6. Carbon 


Pe 


Fig. 6. Effect of partial pressure of hydrogen on carbon 
deposition (“ standard” conditions, except total pressure 
and H, molar feed ratio). 


deposition was found to decrease linearly with 
hydrogen partial pressure, showing the sup- 
pression of coke formation, apparently by 
saturation of olefinic precursors in proportion to 
the hydrogen partial pressure. Conversion, 
however, decreased linearly with increased 


hydrogen partial pressure because of the corres- 
pondingly reduced reactant partial pressure. 
Although the curves shown in Fig. 7 obviously 
apply specifically to the “ standard ” conditions 
and the physical system used in this particular 
investigation, they illustrate very strikingly the 
favourable effect of hydrogen addition on the 
one hand, namely, the diminished coke deposition 
and the associated higher catalyst activity ; and 
the unfavourable effect on the other hand, 
namely, the decreased conversion due to decreased 
reactant partial pressure. 
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Change in conversion vs. average H,/CgH), in 
mixture (all conditions “standard” except 
H, C,H, molar feed ratio). 


Fig. 
reaction 


REGENERATION 


The evolution rates at various levels of original 
total carbon deposition are shown in Fig. 8. As 
will be noted, two distinctly different rates of 
carbon evolution were involved. The primary 
burning occurred at a high rate, and at its 
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Fig. 8. Carbon evolution rates at various levels of original total carbon deposition. 


completion, the rate decreased abruptly, being 
followed by a secondary burning which occurred 
at a much lower rate. This suggested that coke 
was first removed from the more accessible sur- 
faces of the catalyst during the primary burning 
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Fig. 9. Rates of carbon evolution during regeneration vs. 
carbon concentration on catalyst. 


period, and that the remainder of the coke 
burned off at a far lower rate from the less 
accessible surfaces of the catalyst. 

The rate of carbon evolution as a function of 
the instantaneous carbon concentration on the 
catalyst is shown in Fig. 9. The data were 
obtained from a single regeneration four hours 
in duration, which probably implied the removal 
of very nearly all the carbon from the catalyst. 
Thus, values of instantaneous total carbon con- 
centrations were probably reasonably accurate 
representations. Again, two distinctly different 
carbon rates appear to be indicated. 


Primary Burning 

Not much significance could be attached to the 
actual shapes of the primary burning rate curves, 
mainly because neither temperature nor oxygen 
partial pressure was held constant during the 
earlier stages of regeneration. At maximum 
burning rates, catalyst bed temperature would 
increase by more or than 100°F and oxygen partial 
pressure would decrease by about 50%. 

The regeneration time required to achieve 
completion of the primary burning is shown in 
Fig. 10. This plot was obtained from the rate 
vs. time data shown in Fig. 8, the primary 
burning being considered completed at a final 
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carbon evolution rate of 0-025°,, carbon on catalyst 
per minute. Total time to achieve completion 
of the primary burning was directly proportional 
to the total quantity of carbon removed during 
the primary burning. 
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Regeneration time required to complete primary 
burning vs. total carbon deposition. 


Fig. 10. 


The first segment of the plot of carbon evolution 
rate vs. carbon concentration indicates that the 
primary burning had second-order dependency 
on instantaneous carbon concentration, although 
this result was somewhat fortuitous inasmuch as 
catalyst bed temperature was higher during the 
earlier part of the primary burning period, as 
mentioned before. Danr and co-workers [1] have 
indicated such second-order dependency in the 
regeneration of cracking catalysts. 


Secondary Burning 


The segment representing the secondary burning 
suggested first order dependency on instantaneous 
carbon concentration and refers to essentially 
constant temperature and oxygen partial pres- 
sure. The secondary burning rate was also shown 
to have first-order dependency on oxygen partial 
pressure. An approximately fivefold increase in 


carbon evolution rate was obtained by the use 
of pure oxygen rather than air. 

The effect of catalyst activity on regeneration 
beyond the primary burning period was virtually 
negligible. After first regenerating with air for 
the standard 15-min. period, the catalyst was 
further regenerated with pure oxygen for variable 
lengths of time between 0 and 60 min., carbon 
being evolved at a constant rate over this entire 
range. During the 60-min. pure-oxygen regenera- 
tion, 25°, additional carbon was removed, based 
on the carbon removed during the standard 
15-min. air regeneration, yet activity was the 
same as that for the catalyst not regenerated 
beyond this standard air regeneration. This was 
taken as further evidence in support of the 
hypothesis that during secondary burning, carbon 
was removed from relatively inaccessible areas 
within the catalyst pellet. These areas, because 
of their comparative inaccessibility, contributed 
so little to the over-all effectiveness of the catalyst 
that no increase in activity was shown, even early 
in the run. 

In view of these observations, it appeared 
possible that the mechanism of secondary 
burning was diffusion-controlled. The burning 
may have been progressive toward the interior 
of the pellet ; or it may have taken place in the 
smaller and more inaccessible pores. 

Assuming a mass-transfer controlling mechan- 
ism, Seaver [7] has suggested a further analysis 
of this secondary burning by considering the 
evlinder as consisting cither of an infinite slab 
or an infinite cylinder. Burning of the coke in 
the cylindrical pellet was assumed to occur at a 
continuously receding boundary area. Equations 
were derived ; and using these, resistance to mass 
transfer at the surface of the catalyst pellet 
appeared to be negligible, as indicated by a 
correlation with the data. Similarly, all terms 
except those designating the carbon concentration 
on the catalyst were probably constant. Thus 
these equations predicted that the carbon evolu- 
tion rate could be related as follows : 

For the slab, 


dCe 
da Cey 


Ce? 
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and for the infinite evlinder, 
dc 

ie k 1 


Cec 


where 


rate of carbon evolution, per cent carbon 
on catalyst (based on weight of fresh 
catalyst) per minute ; 

original carbon concentration (includes 
primary carbon) per cent, based on 
weight of fresh catalyst ; 

instantaneous carbon concentration, per 
cent, based on weight of fresh catalyst. 


Indeed, plots of dCe dé vs. 1 (Ce, — Ce) and 
dCe dé vs..1 In (Ce, Ce) both resulted in straight 
lines. 
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Fig. 11. Hydrogen concentration in coke vs. degree of 

coke deposition (“ standard” conditions except process 

time). 


Coke Composition 

The composition of the coke deposit as it varied 
with the degree of deposition is shown in Fig. 11. 
The decline of hydrogen content with the degree 
of deposition was anticipated in view of the fact 
that the degree of deposition was varied by 
changes in total processing time, thus not only 
varying the degree of coke deposition, but the 
degree of dehydrogenation of the coke. The 
hydrogen content in the coke was found to 
decrease exponentially with process time. 
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Mass transfer between isobutanol and water in a spray-column 


P. M. Heerrses, W. A. Hoive and H, Tatsma* 
The University of Technology, Laboratory for Chemical Engineering, Delft, Holland 


( Received 11 March 1954) 


Summary--aA study of the rate of mass transfer across a liquid boundary beween a discontinuous 
and a continuous phase, using a binary mixture and a droplet system has been made. 

Because of differences in density and movement of the drops the continuous phase appeared 
to be thoroughly mixed. This resulted in a nearly constant concentration independent of the 
height of the column. 

For the continuous phase a relation for the partial mass transfer coeflicient k, in the form 
k.dg/D = Sh = h has been used. 

Satisfactory correlation of the experimental results could be obtained with this equation. The 
constant A in the above relation appeared to be a function of the relative interfacial velocity 
which in turn is a function of the viscosities of the drops and the surrounding fluids. Following 
the concept of Hapamano this function is taken as »./»,+n4- 

This picture of the mechanism is representative of the situation, when the circulation in the 
drop is not restricted by interfacial tension. 

For transfer to the discontinuous phase the unsteady state equation of Kronic and Brink, 
considering the occurrence of internal droplet circulation, was shown to express fairly adequately 
the experimental results for free rise of isobutanol droplets. Water drops showed a divergence 
from the Kronig-line possibly due to the * floating ” effect of water droplets in isobutanol. 


Résumé— Les auteurs ont étudié la vitesse du transport de matiére, 4 travers une surface limite 
entre deux liquides, vers les phases continue et discontinue en utilisant un mélange binaire et un 
systéme a gouttes. 

Par suite de la différence en densité et du mouvement des gouttes la phase continue se révélait 
intimement mélangée d'ot résultait une concentration presque constante indépendant de la 
hauteur de la colonne. 


Pour la phase continue, on a utilisé une relation, pour le coefficient partiel du transport de 
matiére k,, sous la forme de 
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Cette équation permettait d’obtenir une corrélation satisfaisante avec les résultats expéri- 
mentaux. La constante 4 dans la relation susnommée paraissait étre une fonction de la vitesse 
relative se produisant a la surface limite. Cette vitesse est une fonction des viscosités des gouttes 
et des liquides envirronnants. Basée sur la conception de Hapamarp, cette fonction s'écrit 
comme 


Ne + "a 

Cette présentation du mécanisme est représentative quand la circulation dans la goutte n'est 
pas limitée par la tension a la surface limite. 

Pour la transmission 4 al phase discontinue, on prouve que l'équation de Kronic et Brink 
~ équation qui tient compte de la présence d'une circulation interne dans la goutte — exprime 
de fagon assez efficace les résultats expérimentaux en cas d'un vol libre de gouttes d’isobutanol. 
Les gouttes d'eau présentent une déviation de la courbe de Kronic. Cette déviation est causée 
probablement par le phénoméne de * flottaison * manifesté par les gouttes d'eau dans isobutanol. 


_ ® Part of the thesis for Ph.D. of H. Tatsma, presented at the University of Technology, Delft (1953). 
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INTRODUCTION 


For the study of the mass transfer between two 
liquids a partially miscible binary system was 
chosen, because if equilibrium at the interface 
can be assumed, the interfacial concentration is 
then only a function of the temperature, being 
the saturation concentration. A spray column 
was finally used, because in that case the inter- 
facial area can easily be calculated. Other types 
of apparatus with known areas such as the wetted 
wall column and the horizontal extractor showed 
disadvantages over the spray column. In the 
wetted wall column the film at the wall is easily 
disturbed ; in the horizontal extractor a station- 
ary condition is only very slowly reached in the 
viscous region of flow, moreover small differences 
in density caused layer formation, which hin- 
dered the experiments very considerably. 


The known interfacial concentration makes 
possible the calculation of a partial mass transfer 
coetlicient, the advantages of which are apparent. 


The spray-column resembles technical extraction 
apparatus more closely than the others mentioned. 


The study to be given here can be regarded as 
an introductory study to the problems of real 
extraction. As far as known to the present 
authors, Suerwoop [1] and Licur and Conway 
[2] have used the known interfacial area of a 
spray column in their extraction calculations 
whereas and Wetsn [3] and Lappua 
and Siru [4] have used a binary system. 


Metrnop Apparatus Usep 


In Fig. 1 a diagrammatic sketch of the apparatus used to 
contact the one phase in droplet form countercurrently 
with the other is given. As shown, the apparatus is set 
up for contacting the heavier water droplets falling 
through isobutanol as the continuous phase, The des- 
cription will be based on this type of contact. The appar- 
atus could easily be changed for isobutanol drops rising 
in water as the continuous phase. 


The head tanks for both the water (4) and the isobutanol 
(B) were heavy walled glass containers of 60 litre capacity 
supported 3 metres above the floor. A constant head for 
the flow meters was maintained by the method of Mariotte, 
whereby an open capillary, with its tip placed near the 
bottom of the otherwise sealed container insured that the 


+ To be published elsewhere. 


pressure in the horizontal plane with the tip was constant 
and about equal to that of the atmosphere. 

Provisions were made for filling the vats from supply 
vessels by air pressure and for venting during the time 
of filling. The capillary and vent of the water head tank 
were provided with dust filters and those of the isobutanol 
head tank were connected with a CaCl, drying tube and 
a jar of NaOH in series to prevent the absorption of 
atmospheric moisture in the dry isobutanol. 

By syphoning action the liquids were led via their 
respective capillary flow meters (C and D) and control 
valves to the column distributors. The water stream 
was fed to the upper distribution chamber (2) and then 
to the eight thick walled capillaries (F) which had been 
cut from the same length of capillary tubing. The capillaries 
each 18 centimetres long, with an I.D. of 0-1 em and an 
O.D. of 0-6 em were arranged symmetrically around a 
4 centimetre diameter circle in a rubber stopper which 
snugly fits the experimental column of average diameter 
73cm. A short section of rubber tubing connecting each 
capillary with the distribution chamber, permitted 
individual closure by screw clamps (G). 

After passage through the column, which varied in 
length from less than lem to 128cm, the drops were 
collected by the lower conical section and rolled into the 
narrowed neck where the interfacial position (H) was 
held 1¢em below the constriction by adjustment of the 
overflow leg. The water streamed through this leg and 
by means of a vent at the highest point overflowed 
smoothly into the discharge line. Here a three way cock 
was provided for taking individual samples during the 
collection of the composite sample over a material balance 
period. 

Isobutanol, the continuous phase, was discharged 
directly under the rubber ring joining the lower distributor 
to the column through four evenly spaced headers (J). From 
this annular space, it flowed smoothly and without undue 
disturbance under the circumference and into the column 
proper. 

Provisions for sampling the isobutanol discharge from 
the column were the same as described above for the 
water phase. 

For the investigation of rising isobutanol droplets, the 
whole column was simply inverted and the necessary 
changes made in the flowmeter capillaries for measuring 
the new flows. 

All flow lines were provided with vent valves at the 
highest points to permit easy removal of air bubbles and 
prevent possible blocking of flow. In addition both lines 
contained valves for sampling the inlet streams. An 
Anschutz thermometer, accurate to 0-05°C was placed in 
the line by each flow meter for determining the entering 
stream temperatures and one installed in the column for 
the determination of the temperature under which the 
mass transfer took place. The surrounding room tem- 
perature was kept as constant as possible between 20 
and 21°C, 
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After the start up, a time suflicient for the 
continuous phase to displace four times the 
content of the column was taken to reach dynamic 
equilibrium, as was proved by experience. This 
“run in” time varied from 3-48 hours. During 
this period the head of the flowmeter for isobu- 
tanol and for water the temperature 
of the isobutanol (f;) and water stream (f,) and 
in the column (f,) were read at regular intervals 
and if need be readjusted. After reaching steady 
state conditions the following experimental data 
were taken during the material balance period : 


Fig. 1. 


1. The time of fall of the droplets from each 
capillary (Op in sec) as a mean of 15 determina- 
tions taken from the moment the drop left the 
capillary up to the moment it passed the brim 
of the column. 

2. Number of droplets escaping from a capil- 
lary per second (A,,,). This number was deter- 
mined at low water rates by measuring the time 
necessary to form 20 drops, five times per capillary. 


If A,,,, is larger than 2 a stroboscope was used, 


determining for each capillary three successive 
frequencies at which the drop seems to be at 


Sketch of the extraction apparatus. 


at 
Tk 
1954 
| om 
| 
i 
| | Jk | 
124 


rest, and taking their greatest common measure. 
3. The temperatures ¢,, t,. and the heads 
of the flowmeters AH; and AH, at 2 minutes 
interval (in some experiments 3 minutes). 
4. Five timed volume samples, collected in 
tared flasks, of each exit stream. 


After the end of the experiment, the flow of the capil- 
laries was shut off by closing the screw clamps and then 
immediately closing the continuous phase control valve 
and the overflow leg. The column was then drained at a 
velocity of 2ml sec while at least six representative 
200 ml interval samples were taken of the content. 
Finally a sample of the entering isobutanol stream was 
taken and after reweighing the material balance sample 
flasks, all sample flasks were sealed with paraflin wax to 
await analysis. 


ANALYSIS 


The amount of water in isobutanol was deter- 
mined by means of the Karl-Fischer reagent by 
a slightly modified titration method as described 
by and Surru [5]. The two-solution 
technique was used [6]. 

The amount of isobutanol in water was found 
by determining the carbon content by combustion 
with oxygen over copper oxide at 750°C and 
determining the amount of carbon dioxide formed 
by adsorption in NaOQH-asbestos (Ascarite) and 
weighing. This method will be described else- 
where. 

The water used was distilled water. The 
isobutanol was obtained by distillation of a 
commercial product in a column with sieve-plates; 
boiling point 108-1°C at 760mm Hg, nj = 
1-3959, water content < than 0-05°, 


EVALUATION OF THE DATA OBTAINED 


From the measurements collected the following 
data could be obtained : 


1. The amount of the discontinuous phase 
entering the column (Q,, in em* sec) from the 
mean readings of the capillary flowmeter. A 
correction for deviations of the temperature as 
shown by the mean temperature measured is 
used, 

2. The number of drops per sec (A in number 
sec found as the sum of A.,,, for the eight 
capillaries). 
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3. The initial volume of one drop (V, in em*) 
from V, = ~~ It may be noted that even with 


the same capillary and with the same flow rate 
the drop volume is not always the same. In 
particular for water drops in isobutanol the 
deviations are large. Therefore A has to be 
determined in every experiment. 

4. The amount of the continuous phase 
entering the column (Q,. in cm* sec) from the 
mean readings of the capillary flowmeter. 

5. The amount of the discontinuous phase 
leaving the column (Q,, in cm* sec) from the 
time necesssary to fill the ample flasks of known 
volume as a mean, 


6. The end volume of the drop (V, in em*) 
from Vy = > 


7. The amount of the continuous phase 
leaving the column (Q,, in cm®* sec) from the 
amount in the material balance flask and the 5 
sample flasks, divided by the total time of one 
experiment. 


From these data and with the aid of the 
properties of the mixtures some important 
variables governing the mass transfer can be 
derived. 


(a) Diameter and surface of a drop 

Very shortly after leaving the capillary the drop 
takes the form of an oblate ellipsoid. If « and 8 
denote half of the short axis and half of the long 
axis respectively the surface O of this ellipsoid 
of volume V is given by ; 


By a photographic method the relation between 
V and ~ has been determined both for water 


and isobutanol and is given graphically in Fig. 2. 
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weterdrops in sodutano! 
isobutano! drops in water 


Vinmm? 


Fig. 2. The relation between : and the volume V’ for water and isobutanol drops. 


Assuming a linear relation between Vo and A 


by the method of the least squares the following 
relations have been found : 


z 


These relations are given as lines in Fig. 2. 
With the aid of these equations the relation of 
O to V has been calculated and is given in Fig. 3 
for both * water and isobutanol drops. With 
the aid of these lines the surface of the drop could 
be found for every calculated volume. 

Each drop changes in volume during the rise 
or fall. To describe the mass transport to the 
continuous phase a mean surface O has to be 
used. As will be seen later the conditions for the 
use of Hrxon’s mean surface 


“water” drops 101 — 006438 V 


jsobutanol drops :— = 0-97 — 0-0806 


_ 037? — 


except as far as the constancy of 3 is concerned, 


are fulfilled, although the change in volume is 
33°, at the maximum. The use of an arithmetic 


mean surface O = showed deviations of 


only 0-4°%, at the maximum with the surface 
calculated according to Hixon. For convenience 
therefore the arithmetic mean has been used. 

A mean drop radius (r, in cm) has been intro- 
duced, being defined by O = 477,27. The diameter 
d, = 2r,. 


(6) Total amount of drops, P. 


The total amount of drops present in the column 
is found by adding the products of A... * @.4), 
of the eight capillaries. The total exchanging 
surface in the column is then represented by : 


PO; the mean time of fall or ascent (@) by 
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Fig. 3. The relation of the surface of the ellipsoid drops (O) and their volume (V). 


d= 3 and the mean velocity of the drops 


Li 
- 


vy by vy = ~ in which L is the length of the 


column (em). 


(ce) The column concentration y,. 


During the experiments it could be observed that 
in the continuous phase mixing occurred. At the 
outset, for the calculation of the mass transport 
the concentrations of the entering and leaving 
continuous phase were used to calculate the partial 
mass transfer coefficients. A closer examination, 
however, proved that this was not correct. By 
filling the column for j of its height with iso- 
butanol (with 2°, water) and by bringing in a 
small layer of isobutanol with 0-2°, water coloured 
with Sudan III, which did not mix with the 
denser first layer, and by introducing water drops 


into the column in the coloured layer, within a 
quarter of an hour the two ixobutanol layers were 
thoroughly mixed (same colour intensity). It 
was therefore possible that on this mixing effect, 
no appreciable concentration gradient was 
present, except around the drops. This indeed 
proved to be the case. Therefore, for every 
experiment the column concentration was actually 
measured at different heights, as described. Some 
of the results are graphically represented in 
Fig. 4, for water as the continuous phase (upper 
Fig.) and for isobutanol as the continuous phase 
(lower Fig.) giving the column concentration as 
a function of position in the column. (It may be 
remarked that water is introduced at the top end 
of the column, isobutanol at the lower end, both 
with a concentration = 0). In this figure the 
concentration of the exit stream is also given as 
points on the L =0 abscissa, respectively 
L = 63 abscissa. 
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wn <n Apparently during the “run in” period, owing 

™ - —— to mixing and aided by the way the capillaries 

wougnt | are placed, a cylinder of fairly constant con- 
centration is formed. 

The continuous isobutanol phase in the steady 
state (again aided by the manner of introduction) 
flows along the wall in a thin “* layer,” mixes with 
some of the heavier liquid core, and leaves the 
apparatus at a lower concentration than this 
evlinder-liquid. 

For water, which has a lower viscosity than 
isobutanol, this mixing is more intense, and 
water leaves the column at a concentration 
almost equal to that of the column. 

The concentration of the continuous phase in 
the column (y,) will be given as a mean arithmetic 
concentration calculated from the data measured, 
It will be given in weight percentage of the 
component present in the lowest concentration. 


Fig. 4. The mean column concentration (y,) at different 
heights for water (upper fig.) and for isobutanol (lower 
fig.) as the continuous phase. 


It is clear from Fig. 4 that the column con- 
centration is fairly constant over the whole 
height of the column. For isobutanol as the 
continuous phase this concentration is about 1-25 
times the exit stream concentration and inde- 
pendent of the flow rate (different heights denote 
different flow rates). To this may be added that 
at a low continuous isobutanol rate a horizontal 
plane caused by a difference in refractive index 
could easily be observed in the neighbourhood — oo 984 
of the capillary tips. For water as a continuous weight issbutens! 
phase the column concentration is also constant — pig. 5, The densities of mixtures of isobutanol and water 
and almost equal to the exit concentration, at 25°C. 
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| 


‘ 20 st 


0 25 5 75 10 12.5 
——* weight % water 


(d) The density of mixtures of isobutanol and water 


The knowledge of the densities of the mixtures 
was necessary for the evaluation of the results. 
They were determined at 25-0°C and given 
graphically in Fig. 5 as a function of the amount 
(weight °,) of water present in isobutanol and 
of the amount (weight °,,) of isobutanol present 
in water. 

Circumstances made it necessary for these 
densities to be determined at 25°C. The other 
experiments were carried out at about 20°C. The 
densities found were corrected by means of the 
equation : 
+(1- 
Pmt Pm, 25 Pw. 
in which q = weight fraction of isobutanol, m 
denotes the mixture, w: water, i: isobutanol 
and ¢: the temperature. 

The densities of the pure components at 


Pw, 25 Pit Pi,25 


1.07}— 
waterphase:20 SC 
15 0 25 5 75 10 


Fig. 6. The viscosities of mixtures of isobutanol and water at 20-5°C. 
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weight % isobutanol 


different temperatures have been taken from the 
literature. 


The calculated values are accurate enough to 
be used and a redetermination of the densities 
at 20°C has not been undertaken. 


(e) The viscosity of mixtures of isobutanol and water 


The viscosities of the mixtures of isobutanol and 
water have been determined with an Ostwald 
viscometer and calculated by means of the 


equation: » = ip@— in which » = vis- 


cosity which will be expressed in centipoises, @ 
the time of emptying, p the density, i and 7’ 
constants of the apparatus. 


The constant i’ was calculated from the 


dimension of the viscosimeter ; 7 was determined 
with solutions of cane sugar in distilled water, 
and checked on solutions of glycerol in distilled 
water and on distilled water. 
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Fig. 7. The viscosities of mixtures of isobutanol and water at 25-0°C, 


The viscosities found are presented in Fig. 6 
for 20-5°C and in Fig. 7 for 25°C. 

A linear interpolation between the two tem- 
peratures was used for the temperature of ca. 
20-6°C of the column. 


(f) Mutual solubility of isobutanol and water 


In the calculations it has been assumed that 
equilibrium exists, at the interface between the 
drops and the continuous phase. Therefore the 
mutual solubility of isobutanol and water had 
to be known. The determinations were made in 
two different ways. The water content in satur- 
ated isobutanol solutions was determined by 
shaking isobutanol and water for 24 hours at 
constant temperature and by titrating the 
isobutanol phase after settling. 

For the water phase the temperature of for- 
mation of two phases on heating solutions of 
known composition was determined. The 
criterion was the first formation of a slight 


turbidity. The results obtained are given in 
Fig. 8. For the isobutanol content of water 
slightly lower values than obtained by Micne.s 
[7] have been found. Use of the method described 
for the isobutanol phase gives values lying close 
to those of Micnes. 


(g) The diffusivities of isobutanol and water 

The diffusivities of water in isobutanol (D, in 
sec) and of isobutanol in water (D, in em*/sec) 
were determined with a Schlieren method. The 


apparatus used was that of Wiener [8] as 
improved by Tuoverr [9] and Lamm [10]. 
Assuming that a linear relation exists between 
the refractive index and the concentration, the 
following have been found at 20°C as mean values 
over the whole concentration range involved : 


D,, = (20 + 0-4) 10° cm? ‘sec 
D, = (6-0 + 0-5) 10° em? /sec* 


_ * These values were determined by courtesy of Prof. 
H. Kramers by L. M. Hamaker. 
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Fig. 8. The mutual solubility of isobutanol and water as 
a function of the temperature. 


weight water 


The rather large inaccuracy was caused by the 
fact that the container was not absolutely plane- 
parallel and moreover that the air temperature 
was not thermostatically controlled. 

Measurements with a better apparatus are not 
yet finished. 


(h) Partial mass transfer coefficient in the con- 
tinuous phase (k,) 

The amount of substance (water if isobutanol is 
the continuous phase, isobutanol if water is the 
continuous phase) transported in unit time 
(seconds) from the drops to the continuous phase 
can be found from the concentrations of exit 
and entering streams and the quantities. 

The concentrations are denoted for the con- 
tinuous phase by y in weight %, for the dis- 
continuous phase by 2, also in weight °,. 

The subscript 1 refers to the entering stream, 
the subscript 2 to the leaving stream, subscript ¢ 
to the continuous phase, subscript d to the dis- 


continuous phase, subscript k for the column, 
subscript s to the surface of the drops. Therefore 
tho amount of substance transported in unit 
time, expressed in g /sec is : 


1 
Too (Yo Poe — Ya Pre 


As has been shown the concentration in the 
continuous phase can be taken as constant. 
Denoting the surface concentration by y, (weight 
%,), the driving force for the mass transport in 
the whole continuous phase is : 

1 
100 (y, Ye Pr) in g/em$ 
A partial mass transfer coeflicient for the con- 


tinuous phase therefore in this case can be given 
and defined as : 


Y2 Pre Pic Vic 
PO(y, Pac — Yu Pr) 


in which PO is the surface of all the drops 
present, k, is therefore expressed in em sec. 


(i) Dimensionless groups used 

For the description and correlation of the mass 

transport in the discontinuous phase the dimen- 
sionless groups : 

Ts Psd Pra 


for the continuous phase the Reynolds number : 


Re, = *, the Schmidt number: Se, = 
Ne pe D 
and the Sherwood number: Sh = “Se are 


used. p, and », are determined at the column 
temperature ¢, and the column concentration y,; 
a, and y, at the temperature ¢,. 


RESULTS 


The results of the experiments have been evaluated 
in the manner described. 

To save space the experimental figures will not 
be presented. For readers interested in the actual 
measurements the data are to be had in tabular 
form from the first-mentioned author. 
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Discussion or Resuits 
A. Mass transfer in the discontinuous phase 


To study the results obtained in Fig. 9 for 
isobutanol drops and in Fig. 10 for water drops 
ae is plotted against Pad — Pad 


Vs Psd 71 Prd 
three cases, viz. : 


1. Mass transfer in the sphere by diffusion only 


According e.g. to Newman [11] assuming a 
constant boundary concentration and uniform 
initial concentration the equation of this line, 
which will be called the diffusion line, is : 


— nt nt 


2 


An approximation of this line for a short time 
of contact, can be found by using Hicstr’s 
equation [12] for diffusion from a plane surface 
in an infinite thick layer and integrating over the 
time, which vields : 


6 / De 
€= 
Vr re 


This equation is somewhat easier to handle, but 
gives too high a mass transport because a plane 
area is assumed and the concentration in the 
centre of the sphere is taken constant and equal 


to the initial concentration. For values 


'd 
of 6-10 and 70-10 the differences are 0-1 and 
3°, respectively for the fraction of the possible 
amount taken up. This second line is not shown. 
Moreover the diffusion line is not shown in Fig. 10. 


2. Mass transfer by diffusion and circulating 


Using the stream function of Hapamarp [13] in 
a liquid sphere moving under the influence of 
gravity in a continuous medium, neglecting the 
quadratic velocity terms, Kronic and Brink [14] 
have calculated the mass transport in a liquid 
sphere by diffusion and circulation. The limiting 
assumptions made were : 


1. The sphere has a uniform initial concen- 
tration and the surface concentration remains 
constant ; 


2. The time of circulation (the time which a 


liquid particle following a stream line needs to 
reach its starting point) is small compared with 
the diffusion-time (the time necessary to reduce 


1 
Peg — Poa to (7, — Pra) 


3. As for Hapamarp the interfacial tension 
is zero. 
The equation given is : 
10? 
Pra 
In the original paper only two values for A 
and ,« are given. The first part of the line is 
curved rather strongly, and for this part more 
terms are needed. By courtesy of Prof. R. 
Kronic some more terms have been calculated 
by L. vax peN Branpecer and are presented 
in Table 1. 
Table 1 


1-678 
8-48 
21-10 
38-5 
63-0 
123-8 


For @ — 0, there has been no mass transfer. 
Therefore 
3 @ 
AS = 


As can be calculated : 


> 70-10, 
the fifth term already can be neglected. 
The above equation is plotted in Fig. 9 and 
10 as the Kronig line. 


3. The experimental data 


Comparing the experimental data with the two 
theoretical curves it must be kept in mind that 
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certain conditions postulated in the derivation 
of these lines are not fulfilled in the experiments, 
viz. : 

1. The drops are not spheres ; 


2. The interfacial tension although small (2 
dyne em) is not zero; 

3. The velocity region lies outside the Stokes’ 
region as is shown by Re, and by the ellipsoid 
form of the drops ; 


4. During the formation and during the con- 
fluence mass transfer also takes place (called 
entry-effect and end-effect) ; 


5. The diffusivities used are mean values and 
taken as constant ; 


6. Heat effects accompany the process, 


The results with isobutanol drops —As is clear from 
Fig. 9 with the exception of the points for very 


De 
short columns (low —, values) a curve can 


be drawn through the experimental points. This 
experimental curve lies lower than both the 
diffusion curve and the Kronig curve. This 
means a higher rate of mass transfer caused by 
the entry- and end-effects. An extrapolation to 
De 
ra? 
entry- and end-effects however is not permissible, 
because the experiments were carried out with 
varying amount of drops and varying time of 
drop formation (@,); therefore with varying 
entry effects. 


= 0. in order to determine the sum of 


For 0, + 
> 0, 


loge > — w; for 


log « > 0. 


The end effect will depend among other things 
on the amount of drops. The behaviour discussed 
is confirmed for the lower column lengths. The 
extraction as carried out gives a bundle of lines 
starting at different ordinates (sum of both 
effects) and having a different slope, the line 


axis (small @,) having 
Md 


lying closest to the 
the greatest negative slope. 

This can be seen in Table 2, where the slopes 
of the connecting lines are given as a function 


of @, for each value of the time for two experi- 


ments with different column lengths (L = 3-5 and 
34-8 cm) but with equal @, and equal continuous 
rates. 


Table2 


in sec, slope 


— 28-4 
— 28-6 
— 39-3 
— 


Apparently at higher values of ~ for the 


times of formation used the bundle of curves 
coincide to one curve. 

Although an extrapolation to @ = 0 (L = 0) 
is therefore not permitted, a rough estimate 
(found by extrapolation) indicates that under the 
experimental conditions used for L=0 log 
« = — 0-1, indicating that the sum of entry- and 
end-effects amounts to about 20°, of the maxi- 
mum possible mass transfer. 

In view of the fairly constant entry- and end- 


effects for high values, a decision can be 


Md 

taken as to what is happening in the drops by 
displacement of the experimental line. A vertical 
displacement of the lines from the ordinate 0, 
— 0-1 to 0-0 would include the assumption that 
the drop leaves the capillary with a uniform 
initial concentration and a constant end effect. A 
uniform initial concentration in the drop is 
approximated by short times of formation, causing 
turbulence in the drops. 

A horizontal displacement includes the assump- 
tion that the end effect, which would necessitate 
a vertical displacement, is small compared with 
the entry-effect and that the drop leaves the 
capillary with a concentration distribution equal 
to that in a rising or falling drop. Although this 
is not very likely, it will be clear that without 
further information about entry- and end-effects 
a definite decision about the displacement to be 
undertaken cannot be made. The following 
procedure has therefore been adopted. At the 


highest values of = obtained the slope of 
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Fig. 9. Mass transfer into isobutanol drops. 


the three lines have been determined for equal 
values of ad (vertical displacement) and for 
d 


equal values of log « (horizontal displacement). 
The results are collected in Table 3. 

Including the deviations from theory 1, 2, 3 and 
5 mentioned earlier, and observing that the drops 
are of such a size that the interfacial tension 
cannot influence the movement in the drops very 


strongly (see later), it can be concluded from the 
figures given, that in isobutanol drops ascending 
in a water phase the transport of water in the 
drops takes place by diffusion and circulation in 
agreement with the predictions of Kronic and 
Brink. The greater negative slope of the experi- 
mental line compared with the Kronig line may 
be caused by a different stream pattern caused 
by the ellipsoid form of the drops under the 


Table 3 
Displacement Vertical | Horizontal 
Line | Diff. Kronig Exper. Diff. Kronig Exper. 
| | | | | | . 
=| (80-00 | 230-240 90-100 50-60 
| 
Slope | | | 
— 20 — 60 | | — 20 
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experimental conditions. In the derivation of the 
Kronig line the quadratic velocity terms have been 
neglected. 

The results with water drops—The results with the 
water drops have been investigated in a similar 
way to those with the isobutanol drops. As can 
be seen in Fig. 10 the results with water drops 
seem to be quite different. The results for the 
different column lengths cannot be represented 
by one line. Moreover, as can be seen, for the 
greater column lengths, on increasing the rate 
of the discontinuous phase (and therefore de- 
creasing the time of formation of the drops (@,) 
and increasing the volume V of the drop and 


decreasing the falling time and therefore (= 
d 


the ratio « first increases (as could be expected 
and has been confirmed for the isobutanol drops) 
but then decreases very strongly which is quite 
unexpected. This has been found for both the 
34-1 em as for the 62-8 cm column as shown by 
the drawn lines. Because on decreasing 0, the 
ratio ¢ for both the entry-effect and the exchange 
during the flight increases, the cause of the 
abnormal behaviour must be found in the end- 
effect. During the experiments it has been 
observed that the water drops float for some 
time on the interface before coalescence. For 
the isobutanol drops a similar effect is far less 
pronounced. According to Manasan [15] the 
residence time of falling water drops at the inter- 
face can be 8 sec in a medium of a viscosity of 
4cp. This behaviour of the water drops — in the 
opinion of the present authors caused more by 
the interfacial tension than by viscosity — gives 
the explanation of the effect found. On increasing 
the amount of drops, coalescence cannot take 
place rapidly enough, and the drops are forced 
through the interface with an amount of the 
continuous phase, This amount of the continuous 
phase increases with the amount of the drops and 
causes too high an exit concentration. Therefore 
the mass transfer found will be too high. 

This drag-effect will increase with decreasing 
area of the interface in the neck, and probably 
with increasing »,. On the other hand an inter- 
face of small area has been chosen to decrease 
the mass transport from this interface to both 


phases separated. It may be remarked that very 
great differences in concentration could be 
observed in the water discharge line, while these 
were absent in the case of the isobutanol drops. 

The drag effect decreases also with smaller 
column length, which could be caused by a lower 
isobutanol concentration. Experiments to verify 
this have till now not been very convincing. 

In view of the above it is even dangerous to 
draw an experimental line through the points 
representing experiments with long @, for com- 
parison with the Kronig line. However, for the 
longest column also starting from the largest 0, 
used, in the first instance the ratio « increas*s on 
decreasing @,. Therefore these points have been 
compared with the Kronig line, which for this 
purpose has been displaced in a vertical direction 
and given in Fig. 10 as the dotted line. From this 
it can be concluded that in the falling water drops 
material is also transferred by diffusion and by 
circulation. A very rough extrapolation of the 
experiments last mentioned shows a sum of entry- 
and end-effects of about 35°, of the maximum 
possible, which is 1-8 times that of isobutanol. 
The discontinuous entry-effect—As follows from 
the foregoing paragraph a quantitative knowledge 
about the discontinuous entry- and end-effects is 
necessary for a more detailed evaluation of the 
results. Several attempts have been made to 
measure this effect. 

First of all the behaviour of the liquid in the 
drop during formation was investigated with 
isobutanol drops coloured blue with cobalt 
chloride. Above a certain water content the blue 
colour changes to pale pink. By means of this 
colour change it could be proved that drops 
formed in 1-3 sec or in shorter times showed a 
very marked circulation. If the time of formation 
was larger than 1-5 sec no circulation occurred. A 
marked pink zone was found in the lower part of 
the drop in which a vertical thin blue channel 
lying in the axis of the capillary was visible. The 
upper part was blue with the exception of a very 
thin film at the interface with the continuous 
water phase. During the formation the horizontal 
dividing plane moved upward with a velocity 
dependent on the velocity of formation. 

Moreover on ascending, — blue - coloured 
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Fig. 10. Mass transfer into water drops. 


isobutanol drops gradually changed their colour 
to pink in the whole drop, with a blue ring in 
the equatorial plane which disappeared last of all. 
This is in agreement with the theory of HapamMarD 
which demands that in a sphere with circulation 
a ring in the equatorial plane on a distance of 
= 0-717, from the centre is at rest. This 
again proves circulation in the isobutanol drops. 

Although the time of formation in our experi- 
ments is shorter than 1-5 sec, an attempt to 
calculate the order of magnitude of the mass 
transfer during formation have been carried out 
under the following assumptions : diffusion takes 
place only; the velocity of diffusion is small 
compared with the velocity of growing ; the drop 
is formed by introducing a constant amount of 
liquid (as in the experiments) of constant con- 
centration (¢, in moles cm*) in the centre of the 
drop ; the drop is a sphere with constant boundary 
concentration (¢, in moles cm*). 


Applying Hicute’s equation for the rate of 
transport (R) at the time @ in moles per sec per 


em*®, R = (ce, — ey) = from a flat surface 
7 


with a constant concentration (¢,) into a layer 
of infinite thickness having concentration ¢, at 
time @ = 0, we obtain for the total rate of mass 
transfer (R’) at time@ in moles per sec to the 
inside of a drop : 


o 
R 
On 


As O°? = p, V, = p, @, in which p, and p, are 


constants the integration gives for the rate of 


mass transfer under consideration : 


R = 4p2*(e, — 


The mass transfer to the drop during the time 
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of formation @, therefore is, expressed in moles : 


0 
D 
p,** (c,—¢o) [P 


24 Pr 4 16 D 

This mass transfer can also be expressed as 
(cg — ¢y) Vy. Therefore : 


24 p 


Vp = 


obtains : 


Ts Psa — Pra dy 

In Table 4 the experimental results for very 
short columns for water (1-5) and for isobutanol 
(71-76) have been collected together with the 
values calculated from the above equation. 


C, — Ce 


Table 4 
| | | 
No. of 6,3 in l—e 
eaperi- in sec. ds - 
ment sectem™ found calculated 
7 O24 2-47 0-196 0-0405 
72 O32 #12 0-221 0-0513 
73 0-38 3-52 0228 00579 
74 0-49 3-91 0-236 0-0642 
75 0-89 5-61 0-253 0-0921 
1 O24 238 | 0291 0-0676 
2 O-34 2-92 0-319 0-0829 
3 0-63 416 0-363 0-1182 
4 1-18 5-90 0-380 0-1677 


As in the equation (1 — e) for the experiments 
increases with increasing but not 
linearly. Moreover, the values found are much 
higher, which is easily understood because : the 
end effect still exists; there will be some mass 
transfer in the short column; circulation in the 
drop for the small @, used has been proved ; part 


By expressing p, and p, in d, (diameter of the 
sphere at time @,) and by rearranging one 


of the drop after it has been released adheres to 
the capillary and finally surface oscillations 
caused by the release of the drop may occur. 

A study of the entry- and end-effects must be 
made in a different way. Until now no success 
was met in this respect. Formation of drops 
Imm above the liquid surface and comparing 
those with drops formed 1mm _ beneath the 
surface met with no success at all since a fall 
even over so short a distance in air caused 
deformation of the drops and in all cases a very 
strong turbulence in and around the drops. 
During these experiments in some cases a small 
bubble of air was trapped in the water drop. 
Photographically it was found that the bubble 
was always in the equatorial plane while the 
maximum distance of this bubble from the 
centre was found to be 0-6r,. This proves the 
circulation of the liquid in the water drops. 

Experiments are now being pursued in a 
different direction. 


B. Mass transfer from the drops to the 
continuous phase 


The mass transfer from drops to the continuous 
phase is governed by two factors: (a) diffusion 
from the interface perpendicular to the direction 
of flow and (b) the liquid velocity along the 
surface of the sphere. It may be remarked that 
in our experiments only the immediate surround- 
ings of the drops have to be included in the 
observations because of the mixing effect observed. 

For a moving solid sphere the velocity at the 
interface can be assumed to be zero in respect to 
the sphere. In such a case, in view of the magni- 
tude of the diffusion coefficients, mass transfer 
takes place mainly perpendicular to the surface 
and not over the surface. Therefore the con- 
centration gradient in the direction of flow is 
small compared with that perpendicular to the 
direction of flow. 

For a falling fluid sphere in which circulation 
takes place, the interface moves in respect to the 
centre of the sphere, if this centre is considered 
at rest. This implies that the top of the sphere 
in the direction of flow always meets fresh 
liquid. During the flow past the sphere mass is 
transferred into the stream from the interface by 
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diffusion and a concentration gradient is built up 
in the direction of flow as well as perpendicular 
to the direction of flow. Mass is transported 
mainly over the surface of the sphere in the 
direction of flow and leaves the sphere at the 
bottomside if the velocities are large compared 
with the diffusivities. 

If the circulation in the fluid drop is limited 
by the fact that either », is not small compared 
with », or that the drop diameter becomes smaller 
than the critical diameter d_, defined by Bonn 
and Newton [16] as: 


d., =2 Ji 
(pa — Pe) 


the velocity at the interface in the continuous 
phase will accordingly diminish. In such cases 
the type of mass transfer from a fluid sphere can 
change gradually to that from a solid sphere. 

The partial mass transfer coetlicient for both 
types of mass transfer must depend in different 
ways on Re, and Se,. 

The correlation found so far in the literature 
takes the form Sh = f + h Re) Se*. For a solid 
sphere, f which is the Sherwood number when 
there is no flow past the sphere, can be calculated 
from the theory as being 2-0. 

For a solid sphere Fréssiinc [17] found semi 
empirically h = 0-552, b = }, a = 1.3. For gas 
bubbles Prerers and van Kreve.en [18] found 
f=Oh=11,b= },a=—}. For heat transfer 
from a solid sphere with constant surface tem- 
perature to a continuous medium Rawnz [19] 
derived: Nu = 2 + 0-60 Pr'* which 
using either the Lewis or Chilton-Colburn analogy 
yields a mass transfer equation of the Fréssling 
type. 

For potential flow across a solid (or any other) 
sphere Bousstnes@ [20] derived an equation for 
heat transfer : Nu = 1-13 Pe® 1-13 Pr® 
or an analogy for mass transfer Sh — 1-13 Re®® 
(Pe = Peclet number). 

The same type of equation can be derived by 
using the Higbie concept of mass transfer and 
assuming either potential flow across a_ solid 
sphere or using Hapamarn’s stream function for 
the continuous phase, yielding : Sh = h Re,°* Se°* 
in which A depends on the thickness of layer 


around the sphere in which a concentration 
gradient exists. 

It seems quite safe to neglect f for the regions 
investigated. Moreover, in view of what has been 
said, the power of the Reynolds number can be 
assumed to be 0-5 in first approximation. For 
mass transfer from moving spheres to the con- 
tinuous phase the following equation can then 
be given: Re’ Se". Both h and a 
will depend on the circulation in and therefore 
outside the sphere. A measure for the type of 
flow will be the ratio between the liquid velocity 
at the interface of the sphere in an equatorial 
plane perpendicular to the direction of the 
undisturbed flow, to the velocity of the un- 
disturbed flow, if the centre of the sphere is 
thought to be at rest. 

For the Stokes’ region of flow this ratio can 
be found from the Hadamard stream function 
and appears to be : 


— "e 
+ 

Although the Reynolds number indicates that 
the experiments have been carried out outside 
"Ne 
Ne "a 

governing the magnitudes of both h and a. 

For the experiments to be described however 
only a small variation of the Schmidt number 
was possible in each case. In view of the above 
and of the results discussed later a will be taken 
as 0-5, making / the only variable. An eventual 
variation in a will therefore manifest itself in A. 

In Fig. 11 Sh as calculated has been plotted 
against Re,) Se} with the exception of two 
experiments (see below). Excluding the experi- 
ments in which y, was not measured but only 
calculated, for the three sets of points by means 
of the method of the least squares three straight 
lines are drawn which have the following equa- 
tions : 


isobutanol, L = 62-8 cm: Sh= —137 +-0-98 Re} Se} 


isobutanol, L =34-1lem: Sh= —78 + 0-96 Re} Se} 
water, L=62-8em: Sh=— 44 Re 


the Stokes’ region, the ratio is taken as 


The negative value found for f has no real 
significance, because an extrapolation of the lines 
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Fig. 11. Sh as a function of Re? Sc.4 for the continuous phase. 


found to Re, = 0 is not permissible, because the 
number of experiments and the regions are rather 
small; moreover the line may be curved for 
smaller values of Re, Sc,. 

A better method of determining the constant h 
is to neglect f and to determine h as the slope of 
the line drawn through every point and the origin 
(see below). 

In Fig. 11 some experiments of Guyer and 
Pristrer [21] as carried out with CO, bubbles in 
water have also been included. Their experiments 


enable one to calculate Sh and we = Re, Se, 


if it may be assumed that the resistance in the 
gas phase can be neglected compared with that 
in the liquid phase and moreover that the CO, 
concentration in the bulk of the continuous phase 
is zero. 

In the experiments in question a slight influence 


of y, on Sh can be observed, in this sense that h 
decreases with increasing y, as measured. The 
value for y, is a mean value of the core and of the 
stream along the wall which has a lower con- 
centration. The mean value used for y, therefore 
will be somewhat too low, which will manifest 
itself in a Sherwood number that is too low. 

The influence will be larger if y, is near y,. This 
is clearly seen in the two experiments which have 
been excluded from Fig. 11 (see before). Here 
y, as calculated was 12-9 and 13-8. A fairly small 
increase to 13-7 and 14-1 respectively would bring 
the Sherwood number on the line calculated. 

Also for the continuous phase entry- and an 
end-effect will be present. These effects however 
will be smaller than for the discontinuous phase. 
In view of the concentration distribution in the 
column during formation and during coalescence 
the driving force will not be far from (y, — y,). 


P. M. Heertves et al. : Mass transfer between isobutanol and water in a spray-column 


» 


@ isobutanol 


02 


Fig. 12. The factor 4 as a function of —"*—, 


Near the capillaries a marked concentration 
gradient has been found outside the exchange 
region, moreover the same was the case near the 
inlet of the continuous phase. For isobutanol 
as the continuous phase the stagnant layer 
beneath the inlet had a concentration of 15°,, 
which must be compared with y, = 16-6°%,. The 
product for partial mass transfer and surface will 
be small for the entry- and end-effects compared 
with that for the rising and falling drops. There- 
fore the partial mass transfer as calculated will 
be somewhat larger than that actually existing. 
The relative difference will become larger if the 
column is shorter and will disappear with a 
column of infinite height. 

The equations as found for the Sherwood 
number for isobutanol confirm this. This can 
also be concluded from the experiments using 
water as a continuous phase. Even for the very 
short column this tendency is present, although 
the stability of this phase is very poor owing to 
the influence of the entering and leaving stream. 


7c 


To compare the mass transfer from solid 
spheres, isobutanol and water drops and CO, 
bubbles to a continuous phase in Fig. 12, the 
factor h as described (h = Sh/Re,Sc,4) is 


plotted against —o for the isobutanol and 


d 

water drops and for the CO, bubbles with total 
circulation, as well as the values which would 
follow if in these cases the drops behaved as solid 
spheres in which case h according to the Fréssling 
equation would be 0-552 /Sc,"*. The values used 
in this last case are given in Table 5. 

In Fig. 12 only these experiments have been 


Table 5 


1.2 
| | | 
1.0 - | 
0.8 
| 
| 
06 
| 
| | 
H | | | VOL. 
02 3 
O water 
0 
0 || 0.4 0.6 0 1.0 
| 
Discontinuous | 
phase Se, 0-522 | h 
CO, 570 | 0-19 0-96 
isobutanol 1790 0-16 | O70 
water 23600 | 0-10 | 0-85 
| | 
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Table 6 
Dis- Y A d., d,, 

Lit. continuous Continuous g/cm* in cm in cm 
phase phase | | measured 

16 air waterglass | 107% 0-5 0-33 

16 air golden syrup 91 1-48 0-33 

16 air water 0-018 73 10 0-5 0-33 
21 cO, water 0-015 73 10 0-5 0-35 (?) 
included in which the velocity at the drop surface Ne 


was at the maximum possible, therefore for drops 
with total internal circulation, for which d,> d_,,. 

From the experiments of Guyer and PrisTer 
it can be calculated that an increase in d,, goes 
parallel with an increase in Sh (Re, Sc,)*. Starting 
from 0-49 for drops of 0-18 cm, for drops of about 
0-47 cm and larger Sh becomes constant 
(ca. 0-96). The critical diameter in this case 
would be of the order of 0-35em. Some other 
values as given in the literature are collected in 
Table 6 in which the critical diameters as caleu- 


lated from d= 2 J have also been 
Ap 


collected. 

For the isobutanol drops (0-15-0-18 em) and 
the water drops (0-17-0-22cem) the calculated 
d,,= 02cm. In view of the fact that the 
calculated data seem somewhat high because the 
experiments lie outside the Stokes’ region, com- 
plete internal circulation has also been assumed 
in these drops. 

The results discussed permit the conclusion 
that for mass transfer from CO,, isobutanol and 
water drops to a continuous phase the Sherwood 
number is directly proportional to Re,°°. 

The data given in Table 5 indicate that in 
these cases Sh is also directly proportional to 
Sc.°* and h therefore a real constant independent 
of Re, and Sc,. 

Starting from solid spheres (or spheres which 
behave as such, like a liquid drop in a gas) the 
proportionality constant h based on Sh = h Re? 
Sc.°* increases rapidly on increasing the ratio 

Ne 


Ne "a 
the maximum internal circulation is possible. 


» that is on decreasing », from o to 0, if 


Between zero and a certain value >0 of ; + m 
which might be very close to zero, h will denen 
on the Schmidt number as follows from the fore- 
going. In these cases h = h’ (Se)-™", in which 
h’ and m (varying from 1/6 to 0) are real con- 
stants. Both hand h’ and also m, are a function of 


me For gq = ©, h’ reaches the value 0-552 
Ne "a 
and m the value of 1,6. 


The value of — = at which m becomes zero 


Ne "a 
is not yet known. As said above this value might 


be very close to zero. 

If the condition that the diameter of the drop 
is larger than the critical diameter is not fulfilled, 
the constant h will be lower than that found for 
drops with the maximum internal circulation. The 
limiting case is that the drops become so small 
that no circulation at all will occur, in which 
case the drop will behave as if solid. 


Part of the experimental work has been carried 
out by D. FLemmine, R. E. Lrran, P. L. SLEVEN 
and A. M. A. J. TreLens as partial fulfilment for 
the examination for Chemical Engineer. 

The authors should like to thank at this place 
Prof. H. Kramers and Prof. R. van pe Laar 
Kronic for valuable discussions on parts of the 
work, 

NoraTIon 
a = power of Sc, 
A = number of drops per sec 
A,, = constant in n term of a series 

b = power of Re, 

¢ = concentration in moles /em* 

d = diameter of drop in cm 

D = diffusion coefficient in sec 
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Jf = Sherwood number in stagnant medium 
g = acceleration of gravity, em /sec* 
h, h’ = proportionality constants 
AH = head of flow meter in em of liquid 
i, i’ = constants 
k = partial mass transfer coefficient in cm, sec 
m = power of Sc, 
n = index of refraction 
Nu = Nusselt number 
O = surface of drops in em? 
P = amount of drops 
Pe = Peclet number 
Pr = Prandtl number 
P).P2=constants 
q = weight fraction isobutanol 
Q = amount of continuous and discontinuous phase in 
em? /sec 
= radius of drop in em 
= velocity of mass transport in moles, see em? 
= rate of mass transfer in moles /sec per drop 


d. 
Re, = Reynolds number 
Ne 


| 


Ne 
D 


I! 


= Schmidt number 


Pe 
kd 
Sh — Sherwood number ry 


t = temperature in °C 
= velocity of drop in em /sec 
volume of drop in em* 
x = concentration of substance in smallest amount 
present in discontinuous phase in weight °%, 
y = concentration of substance in smallest amount 
present in continuous phase weight % 
a = half short axis of ellipsoid in em 
8 = half long axis of ellipsoid in em 
y = surface tension in dyne/cm 
7» = viscosity in poise 
0 = time in sec 
#, = constant in n term of a series 
p = density in g/em* 


« = ratio 


Subscripts 
¢ = of the continuous phase 
cap. = per capillary 
er = critical 


d = of the discontinuous phase 
k = of the column 
i = of isobutanol 
m = of the mixture 


8 = at the surface of a drop 

= at temperature 

v = after formation of the drop 
w = of water 

0 = at time @ 

0 = at time zero 

1 = of the entering stream 

2 = of the exit stream 
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Book review 


Fourth Sympcsium on Combustion. Bailli¢re Tindall 
and Cox, London. 1953. xx + 926 pp. 56s. 


Combustion research has been expanding rapidly during 
recent years. Whereas 400 people attended the Third 
Symposium in 1948, more than 700 were attracted to the 
Fourth Symposium at M.I.T. in September 1952. The 
scope of the Fourth Symposium was deliberately limited 
to the physical aspects of combustion, especially flame 
phenomena and turbulence. Despite this restriction, 112 
papers, more than in the Third Symposium, were finally 
accepted and published in this account of the proceedings. 


A new feature of this series of Symposia, is the inclusion 
of nine invited survey papers, which occupy the first 113 
pages of this volume. These reviews increase considerably 


The largest group of papers, 50 in all, falls under the 
heading of laminar combustion and detonation waves. 
Theoretical analysis has become possible for more simple 
flames while there is experimental work on the effects of 
pressure and inert diluents with increased precision in 
the measurement of temperature profiles and burning 
velocities. The papers on premixed turbulent flames show 
how the concept of the wrinkled flame front is being 
refined. In the more applied section, although the main 
impetus comes from problems arising in rocket and jet 
engines, there are about a dozen papers concerned with 
the burning of fuel jets and droplets relevant to furnace 
combustion. 


The amount of space devoted to discussion, about 2%, 
is rather limited; only a fifth of the papers receive any 


the value of the book to the general reader. Their subjects 
illustrate the scope of the Symposium, dealing with limits 
of inflammability ; ignition ; burning velocity measure- 
ments ; laminar flame structure ; instabilities such as 
cellular flames; quenching, flashback, blowoff ; 
turbulent flames ; burning of fuel jets ; and flame stabi- 
lisation by baffles placed in the stream. 


comments but these are usually pertinent and concise. 
The book contains a large number of diagrams and 
photographs. The production and binding are excellent, 
and as the result of a subsidy, the price is exceptionally 
low. It can be warmly recommended to all concerned 


‘ith busti blems. 
with combustion problems J. C. Lee. 


Conference on Oxidation Processes (Amsterdam) 


A Joint Conrerence of the Institution of Chemical Engineers and Chemical 
Engineering Group of the Society of the Chemical Industry in London with the 
Royal Institution of Engineers (Chemical Engineering Group of the Netherlands) 
and the Royal Netherlands Chemical Society (Section for Chemical Technology) on 
Oxidation Processes, was held in Amsterdam on 6th, 7th and 8th May, 1954. The 
papers presented at the Conference are given below. Arrangements for publication 
of these papers, together with the relevant discussions, through the medium of 
Chemical Engineering Science, are now being discussed. 


The BASF-process for production of acetylene by partial oxidation of 
gaseous hydrocarbons by E. BarruoLome, Ammoniaklabatorium der Badischen 
Analin- & Soda-Fabrik A.G., Ludwigshafen a, Rhein, Germany. 


Oxygen processes in the steel industry by W. J. B. Cuarer, The British Oxygen 
Company Ltd., Morden Factory Estate, London. 


The provision of oxygen for industrial processes by A. M. CLark, The British 
Oxygen Company Ltd., Morden Factory Estate, London. 


The development of a process for the manufacture of p-tertiary-butyl- 
benzoic acid by R. M. Coie, A. W. Farrparrn and K. D. Deriine, Shell 
Development Company, Emeryville, California, U.S.A. 


The use of oxygen in gasification by F. J. Denr, West Midlands Gas Board, 
Birmingham, England. 
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The air oxidation of mixtures of benzene and cyclohexane to phenol by 
M. B. Donatp and J. R. Grover, Ramsay Laboratory of Chemical Engineering, 
University College, London. 

Preparation of cumene hydroperoxide by J. P. Forrurs and H. I, WarrerMan, 
Laboratory of Chemical Engineering, Delft Technological University, Holland. 


Nitric acid oxidation design in the manufacture of adipic acid from cyclo- 
hexanol and cyclohexanone by A. F. Linpsay, Imperial Chemical Industries 
Ltd., Dyestuffs Division, Billingham, Co. Durham, England. 


Oxidations carried out by means of vanadium oxide catalysts by J. Mars 
and D. W. van Kreve.en, Staatsmijnen in Limburg, Central Laboratory, Geleen, 
Holland. 


The mechanism of liquid-phase oxidation by G. H. Twice, The Distillers 
Company Ltd., Great Burgh, Epsom, Surrey, F ngland. 


Herstellung synthetsicher Fettsauren durch Oxydation von Paraffinischen 
Kohlenwasserstoffen mit Molekularem Sauerstoff by G. Wierze.. 
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